Weakly- Coupled non-Abelian Anyons in Three Dimensions 
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We introduce a Hamiltonian coupHng Majorana fermion degrees of freedom to a quantum dimer model. 
We argue that, in three dimensions, this model has deconfined quasiparticles supporting Majorana zero modes 
obeying nontrivial statistics. We introduce two effective field theory descriptions of this deconfined phase, in 
which the excitations have Coulomb interactions. A key feature of this system is the existence of topologically 
non-trivial fermionic excitations, called "Hopfions" because, in a suitable continuum limit of the dimer model, 
such excitations correspond to the Hopf map and are related to excitations identified in Ref. [T] We identify 
corresponding topological invariants of the quantum dimer model (with or without fermions) which are present 
even on lattices with trivial topology. The Hopfion energy gap depends upon the phase of the model. We briefly 
comment on the possibility of a phase with a gapped, deconfined Z2 gauge field, as may arise on the stacked 
triangular lattice. 

PACS numbers: 



I. INTRODUCTION 

Recently, it has been shown that three-dimensional systems 
of free fermions can have defects with Majorana zero-modes^. 
These defects display a ghostly remnant of braid statistics: 
even though the defects are free to move in three dimensions, 
there are two inequivalent ways to interchange a pair of de- 
fects. This situation was analyzed further in Ref. 3 , where it 
was shown that the statistics of these particles is governed by 
an extension of the permutation group which was dubbed the 
'ribbon permutation group.' Motions of these defects realize 
a projective representation of the ribbon permutation group 
which endows them with a non-Abelian anyonic character. 
We will call them 3D non-Abelian anyons, although this is 
a slight abuse of the terminology. These results were general- 
ized to arbitrary dimension and symmetry class. 

Refs. •? and 3 both considered systems of free Majo- 
rana fermions coupled to a position-dependent mass term. 
This mass term was treated as a classical degree of freedom, 
with no quantum fluctuations, which begs the question of 
what happens when the mass term is also allowed to fluc- 
tuate. If the mass term fluctuates about an ordered ground 
state, then the defects which carry Majorana zero modes in- 
teract via a linear confining potential. Therefore, it is natu- 
ral to seek a model without long-range interaction between 
defects. In such a model, 3D non-Abelian anyons would be 
the weakly-coupled low-energy quasiparticles of a system - 
truly a higher-dimensional analogue of anyons in the frac- 
tional quantum Hall effect. In this paper, we succeed par- 
tially in this goal by constructing a microscopic model and 
presenting arguments that it has deconfined defects with Ma- 
jorana zero modes. The interaction between defects has the 
power-law decay characteristic of Coulomb interactions. The 
fermionic excitations of this model are the zero modes associ- 
ated with defects, gapped modes associated with excitations of 
bulk fermions, and topologically nontrivial configurations of 
the mass field (the "Hopfions") which carry fermion number. 
In the ordered phase (in which the 3D non-Abelian anyons 
are confined) the Hopfion will be gapped, as in the model of 



Ran, Hosur, and Vishwanath^ (whose terminology we adopt), 
who identified it as a gapped fermionic soliton. We present ar- 
guments, based on two different effective field theory descrip- 
tions, that the Hopfion is gapless in the Coulomb phase. In the 
first effective field theory, the order parameter is scrambled by 
quantum fluctuations so that the order disappears. We rotate 
the fermions to the local direction of the order parameter so 
that, even when the order parameter disorders and the system 
enters a Coulomb phase, the fermionic band remains gapped 
although, we argue, there are gapless fermionic excitations in 
the form of Hopfions. In the second field theory, we gauge the 
order parameter, thereby suppressing its gradient energy. We 
do this in an anomaly-free way by introducing a fourth spatial 
dimension, of finite extent, so that the physical space is one 
surface of the four-dimensional slab. When the order parame- 
ter condenses, the hedgehogs ^^ become magnetic monopoles, 
interacting via a Coulomb force so long as there are gapless 
fermionic excitations - Hopfions - at the other surface of the 
four-dimensional slab. 

In Section|ll| we derive new topological invariants of dimer 
models. The motivation comes from the model of free 
fermions coupled to dimer s studied in Section II of Ref. |3] 
In order to understand the statistics of quasi-particles in this 
model, these invariants will be essential. However, these re- 
sults are explained in a self-contained manner, and they should 
also be relevant to physicists interested in studying more tradi- 
tional dimer models without fermions. While it is well-known 
that these dimer models on a torus have different topological 
sectors corresponding to different winding numbers, we show 
that there are additional invariants present even on lattices 
with trivial topology, such as a cubic lattice with open bound- 
ary conditions. Understanding these additional invariants may 
be important in simulations of dimer models, as these invari- 
ants imply that simulations using plaquette flips will not be 
ergodic even in a given winding number sector. In the discus- 



sion (Section[Vl|), we raise some open problems regarding the 
energy of these different topological sectors which could be 
addressed using quantum Monte Carlo simulations. 



In Section III we add dimer dynamics to the the model of 



free fermions coupled to dimers studied in Section II of Ref . 
|3 (where the dimers were non-dynamical). We argue that this 
model inherits a Coulomb phase from the ordinary (i.e. with- 
out fermions) quantum dimer model, and discuss the statistics 
of the Majorana zero-mode-carrying defects in this phase. In 



Section IV we give two different field theories which, we be- 
lieve, govern the universality class of the Coulomb phase of 
the fermion dimer model. These field theories predict that 
the Hopfion is gapless, and explain how the Coulomb phase 
evades potential obstructions such as anomalies. Finally, in 
the discussion section, we summarize our results and discuss 
open problems. 



II. TOPOLOGICAL FEATURES OF 
THREE-DIMENSIONAL DIMER MODELS 



It is well-known that the Rokhsar-Kivelson dimer Hamil- 
tonian Hrk has different topological sectors on a torus, cor- 
responding to different winding modes of the dimers. How- 
ever, it has additional topological invariants even on a lattice 
with trivial topology, so long as we consider either a finite lat- 
tice or an infinite lattice with the boundary condition that the 
dimers assume a given columnar configuration at infinity. In 
we present a Z2 invariant of the dimer model; 



Section II A 



we refer to a configuration in which this invariant assumes 
a nontrivial value as a "Hopfion" configuration, for reasons 
explained later. This invariant is present in the dimer Hamil- 
tonian Hrk, but its most natural physical interpretation is in 
the coupled fermion-dimer Hamiltonian H described in Sec- 



tion |III| In Section [nCj we show that this Z2 invariant is just 
the parity of an integer invariant with a simple topological in- 
terpretation in the continuum. 

In this paper, we will be discussing dimer models and 
fermion hopping models on lattices, primarily hypercubic lat- 
tices Z^. We will take the lattice constants to be equal to 1 to 
avoid clutter in the formulas which follow. Then, in d dimen- 
sions, we will use bold-faced vectors to denote points in the 
lattice, r G Z^. Sometimes, we will instead use latin indices 
z, j, /c, . . . to denote points in the lattice, assuming some arbi- 
trary ordering of the points. In three dimensions, which is the 
main focus of this paper, we will use x^, i = 1, 2, 3 or x, y, z 
to denote basis vectors of the lattice in the directions of the 
three Cartesian axes. 



A. Hopfion 



To define the Z2 invariant, we introduce an anti- symmetric, 
Hermitian matrix M. This matrix has matrix elements Mij 
with Mij = if i and j are not nearest neighbors, and 
Mij = ±i otherwise. The signs are chosen so that M has tt- 
flux around all plaquettes. That is, if z, j, /c, / are sites around 
a plaquette, then 
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FIG. 1: Illustration of signs in matrix M in two-dimensions. Arrows 
between sites show the signs. An arrow pointing from site i to site j 
implies that Mij = -\-i and Mji = —i. There is tt-Aux around each 
plaquette. 



In Fig. [T] we show an illustration of matrix M for d = 2 
dimensions. Such a matrix M can be found for any planar 
lattice. In order to find this M for a planar lattice, add addi- 
tional bonds to the lattice if necessary to triangulate the lattice; 
then, choose phases for the Mij to put +7r/2 flux in each tri- 
angle, giving each square plaquette tt flux; finally, remove the 
added bonds. On a planar lattice, the phases to do this can be 
chosen inductively, by choosing them on a connected sublat- 
tice which is also triangulated and increasing the size of that 
sublattice by adding one triangle at a time so that it remains 
triangulated. Such an M can also be found for some higher di- 
mensional lattices such as a three dimensional cubic lattice or 
hypercubic lattice in higher dimensions. To find such a matrix 
M on a hypercubic lattice, we proceed inductively: suppose 
we have such a matrix M^^^ on a hypercubic lattice in d di- 
mensions. Then, consider a (i+ 1 -dimensional hypercubic lat- 
tice as stacked hyperplanes, each such hyperplane containing 
a (i-dimensional hypercubic lattice. In each hyperplane, we 
use the matrix M^^\ but we alternate the sign of this matrix 
from one hyperplane to the next. Then, we orient the arrows 
connecting the hyperplanes so that they all point in the same 
direction; i.e., in three dimensions, we point all the vertical 
arrows in the up direction and in each plane we stack matrices 
as shown in Fig. ([T]) with alternating signs. 

Now consider any dimerization pattern on a lattice with no 
defects, so that each site has exactly one dimer leaving it. 
We specify this dimerization pattern by a set of numbers Uij , 
where i and j label sites connected by a bond of the lattice 
as before and Uij = 1 if that bond is occupied by a dimer 
and Uij = otherwise. Then, we define an index uh, taking 



values in {0, 1}: 

n^ = [l-sgn(Pf(7V))]/2 
where the matrix N has entries given by 



iV^jf — iVlijTiij . 



(2) 



(3) 



MijMjkMkiMii 



-1. 



(1) 



The matrix n determines which entries of N are non-zero, 
and M determines whether the non-zero entries are -\-i or —i. 
One can directly check that the sign of the Pfaffian does not 
change under plaquette moves, precisely because of the 7r-flux 
condition ([T]) on the entries of matrix M. In other words, if 
any plaquette containing two dimers is 'flipped', as depicted 
in Fig. [2j then the Pfaffian of N is unchanged. Therefore, uh 
is a Z2 invariant of a dimer configuration. 

As shown by Fortuin and Kasteleyn, the number of dimer 



A crucial question is the energy of a Hopfion. This is dis- 
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FIG. 2: The basic plaquette flip move depicted above does not change 
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FIG. 3: The red configuration of dimers is the simplest dimer ar- 
rangement which has nn — 1. We call such a configuration a "Hop- 
fion", following Ref. Jj The blue configuration has nu — 0. The 
blue pattern is assumed to be repeated outside of this 18 site section 
of the lattice. The arrows are determined according to the the vr-flux 
rule. 



cussed in Sections III II D and whether the Hopfion has a 
nonzero energy or not depends crucially upon whether or not 
the dimers are in an ordered phase. 

As a final remark, it is also possible to define nu on in- 
finite lattices with boundary conditions that the dimers as- 
sume a fixed columnar configuration at infinity. This is neces- 
sary when comparing to topological results on the continuum 
model in the next two subsections. We define the configura- 
tion on an infinite lattice in which all dimers are in a columnar 
configuration to have uh = 0. Then, given any other dimer 
configuration which assumes the given columnar configura- 
tion at infinity, since the two configurations only differ on a 
finite set of sites we can compute uh by the relative sign of 
the Pfaffians on that finite set of sites. 



B. Relation between the Dimer Model and the 0(3) 
Non-Linear a model 

In this subsection and the next, a dimer cover, or equiva- 
lently a dimerization, will mean one without defects. 



1. General Theory 

Let us assign a unit vector to each point r of the cubic lattice 
according to the rule: 



coverings of a planar lattice is equal to Pf(M): every single 
term in the Pfaffian contributes a +1. Stated in our language, 
Fortuin and Kasteleyn showed that nu = for any dimer con- 
figuration on any planar lattice. For a non-planar lattice, how- 
ever, the situation is not as simple, and we can have uh = 1. 
In fact, a bilayer lattice suffices, as the simplest configuration 
with uh = 1 is given in Fig. |3] It has 18 sites, arranged in 
a 3 X 3 X 2 lattice. We refer to configurations such as this 
as "Hopfions", due to their connection (explained later) with 
the Hopf map. The existence of the Hopfion explains why it is 
hard to count dimer coverings of a three dimensional lattice: 
Pf(M) counts the difference between the number of config- 
urations without a Hopfion and those with a Hopfion, rather 
than their sum. 

Note that while uh invariant is indeed invariant under pla- 
quette moves, it is not invariant under permutations of the 
dimers around longer loops. For example, in the pattern in 
the figure, permuting the dimers around the 8 sites along the 
outside of one of the squares in a given layer changes h. Thus, 
one may wonder how relevant this invariant is for the physics 
of the dimer model. After all, in any experimental realiza- 
tion of a quantum dimer model, there will likely be some am- 
plitude for longer loop moves. Also, in the original motiva- 
tion for the dimer model as an approximation to the behavior 
of spin- 1/2 systems, there also is some amplitude for longer 
loops moves. However, in the coupled dimer- fermion model 
which we will introduce in the next section, we will see that 
the invariance of uh is protected by superselection rules. 
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The vector points either in the direction of the dimer which 
touches that site or the opposite direction, with the sign 
alternating from site to site in the direction of the dimer. 
In this way, we can represent a dimer configuration by a 
unit vector field which points along one of the axes or, in 
other words, by a map from the lattice to the octahedron 
niat : ^^ -^ Octahedron. More generally in d dimensions, 
if the dimer lies parallel to the z-th direction, 1 < i < (i, as- 
sign the vector ± ^ at the two ends of the dimer according to 
the rule: + if the z-th coordinate is even and — if odd, giving a 
map Z^ — > d-Octahedron. Since the two assignments agree 
at both ends of the dimer, we extend the vector field ± -^ to 
be constant on the dimer. To connect to the theory of topologi- 
cal defects, we would like to go further and define a canonical 
smooth extension ncont • ^ — > S^~^. 

The potential difficulties in constructing such an extension 
are (1) that there may be multiple possible extensions locally 
so that the map from niat to ncont is a one to many map (i.e. 
the extension fails to be canonical) and (2) that some configu- 
rations niat may necessarily have singularities in their exten- 
sions to ncont. i-G- points where we must have ncont = 0- 
Rather than writing an obscure formula for ncont, in the fol- 
lowing paragraphs we discuss the issues encountered in con- 
structing ncont and their resolutions. For now, we set (i = 3 
and make the relevant extensions shortly. We view the octa- 
hedron as a discrete approximation for S'^ . Let us imagine 
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FIG. 4: Examples of dimer configurations on a plaquette and the 
corresponding simplices of the octahedron. 



inscribing an octahedron inside a sphere, as shown in Figure 
[5] and then radially projecting it onto the sphere. Now con- 
sider any plaquette in the lattice. The four corners of this pla- 
quette are mapped to four (not necessarily distinct) vertices 
of the octahedron. So long as niat at neighboring points does 
not point in antipodal directions, call such a dimerization cau- 
tious. Moreover, we call a dimer covering very cautious if 
Xi and —Xi dimers never touch sites lying on a minimi lat- 
tice (i-cube. A useful way to describe this is to say that the 
coordinate- wise L^ distance of points p and q on the lattice 
must be greater than one if p and q touch oppositely-directed 
dimers. If the dimers are very cautious, then the piecewise 
linear chain connecting the four corners of a plaquette unam- 
biguously bound a simplex (point, edge, or face) of the octa- 
hedron. Some examples are given in Figure |4^,b,c. 

We now assume our dimer covering b is very cautious and, 
since the construction is general, we work in (i-dimensions 
with U^, S^~^, the d-cube, and the dual d-octahedron, both 
of which are inscribed in S^~^. We want to avoid all arbi- 
trary choices so that our constructions will work in param- 
eter families of dimer coverings. The simplest construction 
seems to be a two step process. First, again using the L^ 
distance, "thicken" each dimer into a rectangular solid. Ex- 




FIG. 5: The octahedron is viewed as a discrete approximation to the 
sphere. 



pand each of the rectangular solids in all directions until R^ is 
fully packed by rectangular solids. Let f{d) be the step func- 
tion (multi-valued at interfaces) which takes each rectangular 
solid to the vertex of the (i-octahedron to which its core dimer 
has been previously assigned. Note that S being very cautious 
implies the /-image of each fundamental d-cube of Z^ lies in 
a simplex of the d-octahedron. Second, smoothen out f{S) by 
convolving it with a smooth function of very small support, 
e > 0, to produce g. The convolution uses the convex struc- 
ture of each simplex of the (i-octahedron. This defines a unit 
vector field ncont on R^ corresponding to any very cautious 
dimer configuration. 

As an aside, we have checked that for d = 3 it is sufficient 
merely to assume S is cautious in order to construct a canoni- 
cal extension, however the construction is somewhat different. 
In particular, apply the previous paragraph only to the plaque- 
ttes in R^ to get an extension with numerous cubical holes. 
The boundary of each hole has six faces, and a combinatorial 
argument shows the image of such a boundary under / does 
not cover the octahedron (6 < 8) and in fact is a contractible 
subset Z of the octahedron. Compressing the image of the 
boundary of a hole under / radially away from c — centroid 
{S'^\Z)ioc = antipote(c) defines a canonical extension over 
the hole. 

The second step is a refinement trick which sends the gen- 
eral dimer cover to a very cautious dimer cover. We have al- 
ready constructed ncont for very cautious dimer coverings and 
indeed a problem appears to arise when n points in antipodal 
directions at neighboring points: there are many possible ways 
to interpolate between n and — n, and it is not clear which one 
to choose. Such a situation occurs when the dimers are in a 
"staggered" configuration, such as depicted in Figure |4] (d). 
Our solution is to imagine that there there is a more refined 
lattice with 1/3 the lattice spacing (in any dimension) of the 
original lattice so that the physical lattice is a subset of the 




FIG. 6: An example of the refinement of the dimer configuration of 
a cube. A configuration which is dimerized in the \xi\ direction is 
inserted for the central cube. 



refined lattice. For notational simplicity we temporarily re- 
vert to 3-dimensions, although refinement clearly applies for 
all dimensions d. If the points of the refined lattice are in- 
dexed by three integers {m^n^p), then the points of the physi- 
cal lattice are the points (3r, 35, 3t), where m, n, p, r, s, t G Z. 
We now make the following assignment of dimer s on the re- 
fined lattice. If the original lattice has a dimer between sites 
(3r, 3s, 3t) and (3r + 3, 3s, 3t), then the refined lattice has 
a dimer between sites (3r, 3s, 3t) and (3r + 1, 3s, 3t) and a 
dimer between sites (3r + 2, 3s, 3t) and (3r + 3, 3s, 3t). The 
dimer on the original lattice has become two dimers on the re- 
fined lattice. If there is no dimer between sites (3r, 3s, 3t) and 
(3r + 3, 3s, 3t) on the original lattice, then the refined lattice 
has a dimer between sites (3r + 1, 3s, 3t) and (3r + 2, 3s, 3t). 
The analogous rule holds for links lying along the y and z 
directions. 

So far our refinement procedure takes a complete dimer 
covering into one that is incomplete. This can be rectified by 
dimerizing each central plaquette in any way which does not 
reintroduce an incautious pair, i.e. antipodal vectors on sites 
at distL°o = 1. This is always possible since each plaquette 
of the unrefined lattice has at least two unfilled opposite sides 
which define the permissible directions for the dimers on the 
central plaquette in the refined lattice. For concreteness, we 
order the coordinates \x\ < \y\ < \z\ (or, more generally, 
\xi\ < \x2\ < \xs\ < '■' < |xrf|) and dimerize central pla- 
quettes in the lowest possible direction. Next, central cubes 
(and then central 4-cubes,. . . , (i-cubes) can all be dimerized 
in the \xi\ directions (see Figure [6]). This completes the con- 
struction of dimer refinement. 

Notice the general dimer cover, after one step of refinement, 
becomes very cautious. Thus for general dimer coverings, the 
composition ncont o (refinement) defines a mapping, call it Vd'- 



The right hand side of the above mapping is a topologi- 
cal space with the C^ compact-open topology, while the left 
hand side is a discrete set. To better compare the two, we 
would like to endow the set of dimer coverings of I/^ with a 
topology. It is customary to think of the set of dimer config- 
urations as a graph, with dimer configurations at the vertices 
and links connecting dimer configurations which can be con- 
nected by the plaquette flip depicted in Fig. [2] However, we 
can go further and promote it to a cell complex Vd by at- 
taching j-cells, j = 1, 2, 3, . . . , to the discrete set of dimer 
configurations. There will be primitive cells: 1-simplices (the 
links of the graph mentioned above), 2-simplices, 3-simplices, 
. . . , d — 1-simplices, corresponding to: plaquette flips, cube 
flips, . . . , (i-cube flips, and further cells corresponding to arbi- 
trary products of all such flips wherever these flips are realized 
disjointly. A cube flip is composed of two plaquette flips, as 
depicted in Fig. [7] A 2-simplex is a triangle which we asso- 
ciate to a set of three cube flips, such as the set depicted in 
the diagram on the right of Fig. |7] Note that the sides of this 
triangle are not attached to three links in the graph. This is 
because each arrow in the cube flip is a composition of two 
plaquette flips on an opposing pair of faces. Therefore, the 
corresponding 2-simplex should, in order to preserve symme- 
try, be attached to the diagonal of the square corresponding to 
the product (i.e. disjoint occurrence of) those two plaquette 
flips (i.e. to a link which is not present in the graph because 
it corresponds to two plaquette flips), as shown in Figure [8] 
Higher simplices are defined analogously. 

As one moves across these primitive cells and their Carte- 
sian products, starting with the map to A^^^ on the boundary, 
one inductively interpolates (within Md) the extension of rd 
across each cell. The result is a continuous map 



rd 



Vd 



M, 



Now we restrict to dimer coverings V^ which are colum- 
nar in the +x-direction near infinity (i.e., +xi -columnar ex- 
cept for finitely many dimers) and simultaneously to M^, the 
space of maps R^ — > S^~^ which takes a neighborhood of 
infinity to (1, 0, . . . , 0) G S^-^. This yields 

r^ P^ ^ Ml 

If the 0(3) non-linear a model is to capture the physics of 
the cubic lattice quantum dimer model, then the topology of 
the space of dimer coverings which are columnar in the -\-x 
direction near infinity must be the same as the topology of 



plaquet flip 



cube flip 



Td : {dimer coverings of Z'^} 



>Maps(R'^ -^ 5^-^) 
:=M{d). 



FIG. 7: A plaquette flip is depicted on the left. A cyclic sequence 
of cube flips is depicted on the right. Each arrow corresponds to a 
single cube flip. 




FIG. 8: The 2-simplex is associated to "cube flip." 



unit vector fields which are equal to x near infinity. A more 
precise way of expressing this is that we need rg to be a ho- 
motopy equivalence. More generally, one can ask whether r^ 
is a homotopy equivalence for all d > 1. 

We do not know if this is true. We can, however, prove two 
theorems, both of which are physically significant. 

Theorem 1. r2 : ^2 — ^ -^2 ^^ ^ homotopy equivalence. In 
fact, both D2 and M2 are contractible. 



,0 



Theorem 2. r^ 



inverse s 







T^'a 



K 



M^^ has a weak right homotopy 
^d • ^^^d — ^ ^d' ^•^- ^^ o 5^ c^ id^Q for all d > 1. 
Weak means that technically mapping properties of the target 
space TW^ are only tested over finite-dimensional complexes 
mapped into M^^. 

Corollary. For all k > and d > 1, the map 

(r2)/e : 7r/c(P2) — > 7r/e(A1^) is onto. Thus, the dimer 
space has ''at least as much topology " as the mapping space. 

Proof. Since 7r/e is a homotopy functor, (5^ o r^)/. = id^^po), 
so (rj])/e must be an epimorphism. D 

It seems possible that the r^ are actually homotopy equiva- 
lences (which would imply 7r/e(Pj]) ::^ Tik{M^^) for all k^d> 
1) but we could not prove this. However, a modest extension 
of the proof we will give for Theorem 2 shows that if we let 
V^ be the direct limit of the refinement sequence 

,j^Q refinement ,^.0 refinement ,^.0 refinement 

^d ^ ^d ^ ^d ^ • • • 5 

then in fact r^ extends to a (weak) homotopy equivalence: 

Theorem 1 will be proven in the following subsection on 
d = 2 dimer coverings. For completeness, note the following 
trivial analog of Theorem I in d = 1: V^ and M^ are each 
single points with rj the only possible map. 

We remark that M^ is homotopy equivalent to 
Maps(5^52) and 7^fe(Maps(5^52)) ^ 7Tk^s{S^). The 

homotopy groups of S'^ are well studied and completely 
computed up to 7r64(5'^). This wealth of information trans- 
lates directly via Theorem 2 to detect /c-parameter families of 
dimer coverings of Z^ (columnar near oo) modulo (simulta- 



neous) plaquette flips for all k, such that < k < 61. We 
hope there are creative ways to use this wealth of information 
in condensed matter. 

Let us now sketch a proof of Theorem 2. For a dimer 
covering 6 of Z^, the continuous map m : R^ — > S^~^ is 
characterized (up to irrelevant choices from a contractible 
space) by the hypersurfaces - we call them walls - marking 
the frontiers in R^ between the "colors" -\-xi, —xi, -\-X2, 
— X2,. . . , -\-Xd, —Xd- These colored regions of R^ are the 
preimages of the {d— 1) -dimensional faces of the (i-cube dual 
to the vertices of the c^-octahedron also regarded as projected 

Definition. Let M^^ C A1^ be the subspace of maps 
corresponding to domain wall configurations so that points p 
and q with opposite colors Xi and —Xi (i.e. m{p) and m{q) 
lie on antipodal faces) must satisfy dist^oo (p, g) > t. 

Lemma. For all t, the inclusion M^^ C TW^ is a homotopy 
equivalence. 

Proof. Given m : R^ — > S^~^, there is a radial expansion 
t{m) : R^ — > R^, t{m){p^ 6) = {i{m)p^ 0), for an appropri- 
ated monotonically increasing function i{m) : R+ U {0} 
R+ U {0}, so that the compositions m o t{m) G M)^ 
Since i{m) can be chosen continuously in m, this defines 



^o,t 



o,t 



the required (weak) homotopy inverse, i : M^ — > M^ 
i{m) = m o t{m). D 

The (weak) right homotopy inverse will be a composition: 
5^ = 5^ o i, for a map s^^ which we construct next. As we 
define s'^ the reader may notice a curious ambiguity which 
relates in a precise way to the cell structure on V^. 

Each bond of Z^ has a type equal to ±Xi, determined by 
the direction (z) in which it lies and the sign: + (— ) if its 
smallest i coordinate is even (odd). Each plaquette has a type 
consisting of the set containing the two distinct types of bonds 
in its boundary. Similarly, for all unit /c-cubes of Z^, 1 < A: < 
d, culminating with a set of d colors representing the type of 
a d-cube. 

Pick t > 2 and m' = i{m) G A1^'^ Let us begin by 
defining an over-complete dimer covering 5{m'). For each 
closed unit lattice, /c-cube c of Z'^ for 1 < k < d, we put 
dimers on all the bonds of c precisely when the set of col- 
ors (i.e. m'~^ (faces of S^~^)) present in c is identical to its 
type. Among such /c-cubes, we call those which are maxi- 
mal under inclusion active. To make S{m') merely complete, 
polarize each active /c-cube in the direction of its centroid's 
color. There are two things to notice about this rule. First, the 
rule never produces conflicting instructions from two active 
/c-cubes, k = /ci, /c2, for the reason that two distinct active 
/c-cubes, k = /ci, /c2, never intersect. To see this, notice that 
if two /c-cubes, k = /ci, /c2, intersect, the requirement that op- 
posite colors have distance > 2 implies that the union of the 
two color sets is also a consistent color set (i.e. no antipo- 
dal pairs of colors). Thus the span of the two /c-cubes would 
itself have been a larger active /c-cube, thereby contradicting 
maximality. The second thing to notice is that the rule is not 
continuous and does not always define a unique dimerization. 



The two are actually aspects of the same phenomenon: as a 
domain wall between colors sweeps across the centroid c of a 
/c-cube, there will be instances where c acquires two or more 
colors. It is precisely the role of the primary cells of V^ to 
provide the target space to interpolate between the dimeriza- 
tions defined by the potential polarizations of an active /c-cube 
when c lies in multiple color domains. The product cells pro- 
vide a similar target space for disjoint, simultaneous crossing 
of domain boundaries by two or more active centroids. Tech- 
nically the above discussion only defines s^^ on the strata of 
finite codimension. This is why we are only constructing a 
weak right homotopy inverse. This completes the description 
of5',:X°/^D0fort>2. 



To prove Theorem 2, it suffices to show that r^ o s'^ot 



rSofiO 



i&j^o , i.e. that the composition is homotopic to 
the identity. It is sufficient to show for m' G V^^ that rj] o 
s'^{m'{r)) and m'{r) are never antipodal, r G M^, or even 
oppositely colored in S^~^, for then the two are joined by a 
canonical homotopy. But, consequent on our definitions of r^ 
and s'^, r^ o s'^{m' {f)) will have the color of m'{f) for some 
r where dist^oo {f^r') < 1. 

Since t > 2 > 1, the two maps are never antipodal, com- 
pleting the proof of Theorem 2. 

In the direct limit or stable context ( ^ ) the preceeding ar- 
guments localize and the opposite composition also becomes 
homotopic to the identity. For p-\-q = d,p and q non-negative 
integers, let P^^ be the space of stable dimers on Z^ periodic 
in g'-coordinates and +Xi -columnar near oo (unless q = d), 
T^ the g-torus, and Maps((T^ x Bp, T^ x S'^-^); (5'^-\ *)) 



are denoted M^^. 



Theorem 3. There is a weak homotopy equivalence Vp 
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2. Special Knowledge When d = 2 

Although our emphasis in this paper is on d = 3, we have 
developed the general theory, and so will also explain how 
the Conway-Thurston"^ "height function" allows a parallel ap- 
proach for d = 2, which can yield additional information. 

Dimer coverings are, of course, the same as domino tilings. 
Domino tilings may be understood as follows. For any group 
G and set of generators S, we can define the Cayler graph 
which has a vertex for each element of G and a link between 
any two vertices g, g' if g = g' s for some s ^ S (with dif- 
ferent colored links corresponding to different elements in S). 
For instance, the square lattice is the free group F{x^y) on 
two letters x, y modulo the relation [x^y] = xyx~^y~^ = e 
(e is the identity element). We can identify a square plaque- 
tte at the origin with the relation xyx~^y~^ since it contains 
links corresponding to x, y, x~^, and y~^. We can iden- 
tify any other square plaquette with a conjugate of the re- 
lation gxyx~^y~^ g~^: from the origin, it goes first to ver- 
tex g, follows the links corresponding to x, y, x~^, and y~^, 
and then returns to the origin. Now consider F(x, y) with 
the two relations ri = [x'^^y] = e, or r2 = [y'^jx] = e. 
The quotient G is a "Conway group". The Cay ley graph of 



G is a three-dimensional graph. Now consider a domino at 
the origin which lies in the x-direction. We can identify it 
with x'^yx~'^y~^ since it is twice as long in the x-direction 
as the y direction. Similarly, a domino at the origin which 
lies in the y direction can be identified with the other relation 
y'^xy~'^x~^. Furthermore, any domino can be identified with 
a conjugate g x'^yx~'^y~^ g~^ or g y'^xy~'^x~^ g~^ of one of 
the relations. In this way, any domino can be lifted into the 
three-dimensional Cay ley graph of G. Any domino tiling lifts 
to a "rough surface" in the graph. 

The group G fits into the following (non-split) short exact 
sequence: 



h 



^D 



z : 


^^G ^ 


^ z^z- 


-^ 


UJ 


l\\ 


l\\ 




Tc^y] 


{xM[x\yUy\x]} 


{x,y\[x,y]} 





This short exact sequence can be understood as follows. The 
first group, Z, gives the possible values of the height function. 
The second group, G, enapsulates the domino tilings. Each 
element of G corresponds to a particular domino in a domino 
tiling of the square lattice. The map h gives the height of 
that domino. The third group, Z Z is the square lattice. 
There is no natural map back from Z^ = Z Z back into 
G, however if the plane is tiled by horizontal and/or vertical 
dominoes in pattern D (and the origin is marked), then this in- 
formation defines a set-theoretic splitting £d (not a homomor- 
phism) as above. The values of i^ mod 4Z are independent 
of D, but £d itself reflects the tiling. Similarly, there is a set- 
theoretic splitting h (above) determined by the condition that 
7~^[a, 6] 7 — > ±4, the sign depending on whether the total 
number of symbols a, b, a~^, b~^ in 7 is even or odd. The 
composition h o i^ is the "height" function on D. It can be 
described as follows. Checkerboard color the 1 x 1- squares of 
the plane. Starting at e = (0, 0), move (in any way) along the 
edges of the dominoes to a site p ^ I? . There are four pos- 
sibilities as you travel along each length one lattice bond and 
for each one add a term ±1 according to the rules: traversing 
a bond with a white (black) square on your left add +1 (—1). 
The sum of those signs is the height at p, given D. The first 
observation is that a dimer flip D — > D' acts on the graph of 
the height function by laying on, or cutting away, a certain ZD 
"body" from the graph. In the case of dimers on the honey- 
comb, the body is a cube and these pictures are well-known. 
For the square lattice, the body is also convex and is pictured 
in Figure |9] 

The second observation^ (see caption to Figure [s]) is that 
with a fixed loop 7 in I? as boundary condition, interior 
dimer coverings (if they exist) correspond to a discrete 
version of the Lipschitz functions (Manhattan, or L^ -metric) 
with Lipschitz constant < 1. In the case of dimer covers 
fixed to +x-columnar near infinity this becomes: "Lipschitz 
functions consntant near infinity with Lipschitz constant 
< 1." In Ref. 4 an algorithm for constructing the unique 
"lowest" dimer cover filling of a boundary condition is given 
(provided a dimer cover exists), and it is shown that every 
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FIG. 10: The height function of an xi columnar configuration of 
dimers. 



FIG. 9: (Color online) A 3-dimensional rendering of two possible 
height functions defining the top and bottom half of the body. 



3. Periodic Boundary Conditions 



dimer cover is connected to the lowest one by a sequence of 
flips. In fact, the most direct (monotone) sequence of flips 
amounts to removing "bodies" until the graph is lowest. This 
process is canonical except for the occurrence of disjoint 
(hence commuting) moves where two or more bodies are si- 
multaneously removable. In the language of this paper, there 
is a dimer space V2 associated to filling 7 with dominoes that 
has vertices for each such filling, edges for flips, and /c-cubes 
for k simultaneously flips. The arguments of Ref. 4 actually 
show: 



Theorem A. For all 7 either V2 is empty or contractible. D 

In the non-compact case, where the sharp boundary 7 is 
replaced with xi -columnar near 00, dimer covers correspond 
to 1-Lipschitz functions approximately constant near infinity 



as seen from the height function in Figure 10 



While there is now no "lowest function," any compact fam- 
ily K of dimer covers, each a finite alteration of xi -columnar, 
can be canonically lowered to a lowest function constant 
outside a sufficiently large compact set 5 C M^, where S 
depends on K. Thus we have: 



Theorem B. P9 



*, that is V2 is contractible. 



D 



Note that, in a coarse sense, the relation between the 
Conway-Thurston height function h and our vector field v := 

rd{S) is 'U = (sin 27r/i/4, cos 27r/i/4). 

Since the mapping space M2 = Maps((6'^, *), {S^^ *)) is 
also contractible. Theorem B implies Theorem 1 of the previ- 
ous subsection. D 



The height function method can be adapted to periodic 
boundary conditions (torus) and partially periodic boundary 
conditions (cylinder). Let X be either a torus or cylinder. The 
new fact in this context is that now h will be multi-valued 
along essential closed trajectories, h is now a twisted func- 
tion realizing a plaquette-flip invariant holonomy representa- 
tion p : TTi (X) — > R. In other language h is now a section of 
a trivial real line bundle equipped with a possibly non-trivial 
flat connection. 

Picking a base point vertex xq lying on the boundary of a 
domino, xq e X, there will be a unique lowest (relative to 
the connection A) h with h{xo) = 0. This requires check- 
ing that the lowest tiling in the neighborhood of xq is defined 
by the unique model shown in Figure 4.4 in Ref. 4 If X 
is non-compact, a fixed neighborhood of infinity can be speci- 
fied to be xi -columnar and then the notion of lowest is defined 
with respect to this restriction. Arguing similarly to the planar 
case with boundary condition 7, one may see that each real- 
ized representation p corresponds to a connected component 
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of the dimer space V{X) naturally associated to X. Figure 
computes V{X) for X = 4 points arranged as a 2 x 2 square 
torus. There are eight dimer coverings organized into: one cir- 
cle containing 4 (with slope=(0, 0)) and 4 point components 
each with different non-trivial slopes. 

We expect for general tori (and cylinders) X the result will 
be analogous. Components of P(X) of maximal slope will 
be contractible while the other components will be homotopy 
circles (corresponding revolving the average angular direction 
of ncont around the target circle). This will certainly be true in 
some coarse sense but may be true on the nose. A promising 
tool to study this question would be a combinatorial Lapla- 
cian Ah of h. One would like to study dimer coverings by 
deforming to the harmonic representatives, but we leave this 
to later work. 



Slope: 
(0,4) 




Slope 



FIG. 1 1 : The octahedron is viewed as a discrete approximation to the 
sphere. 



C. Continuum Description of Hopfions 

In the previous subsection, we have argued that dimer 
model configurations can be viewed as configurations of a 
field whose target space is the octahedron. If the octahedron, 
in turn, is approximated by the sphere, then the dimer model 
can be viewed as a lattice regularization of the 0(3) non-linear 
<j model. We now make a natural relation between the invari- 
ant uh and the so-called Hopf invariant of the 0(3) model. 

Configurations of the 0(3) non-linear a model are defined 
by maps n{x) from three-dimensional Euclidean space R^ 
to the sphere S'^. Let us consider maps which are constant 
at infinity, e.g. let us consider field configurations such that 
n{x) ^ z 8is\x\ ^ oo corresponding to the boundary condi- 
tion on dimers that they assume a fixed columnar configura- 
tion at infinity. Then we can add a point at oo so that we have 
a map from S^ to S'^. The homotopy classes of such maps are 
in one-to-one correspondence with the integers, 7r3(5'^) = Z. 
The integer A^Hopf which indexes these homotopy classes is 
the linking number between the pre-images of any two points 
on S'^ (the pre-images are closed curves) which may be com- 
puted from the integral formula: 



N^ 



H 



4^7 



d^X CauX da ft • dufl X dxfl 



(5) 



where a^ = ^{d^z 



z^GZ. Alternatively, we can define a^ by 



dyaa - dya^ = -— n • dyU x d\n\ 

47r 

the ambiguity in defining a^ from this relation does not af- 
fect the integral. A configuration of a unit vector field n with 



Nh = 1 is depicted in Figure 12 




z — z^dyz) if we write ft in the form 



The resemblance between the configurations in Figs. [3] and 
12 and the fact that the configuration of Fig. [3] maps to a vec- 
tor field with uh = 1, are not accidental: 
Theorem 4. 

uh = NHmod2, (6) 

Proof. We have shown that different dimer configurations 



FIG. 12: A unit vector field n with Nh = 1. 



with different values of Nh cannot be connected by dimer 
moves consisting of plaquette flips. We have not shown that 
Nh is the only obstruction to connecting two different dimer 
configuration. However, theorem 3, for (p^q) = (<i, 0), says 
precisely this when "dimer refinement" is added to "plaquet 

flip" as an allowable move. In particular, ttq (^3,0) — ^ is 
generated by the dimerization S in Fig. [3] For a general space 
TTo is merely a set, but ttq (^3,0) i^ ^^ abelian group. The 
"sum" of two configurations is represented by taking two balls 
containing the nontrivial, i.e., non-xi -columnar part of each 
configuration respectively and shifting the two balls relative 
to each other by some large even vector V = (2xo, 2?/o, 2zo) 
so that they become disjoint, and finally completing the dimer 
covering outside the two balls by xi -columnar dimers. This 
definition mimics the abelian operation on 7r3(6'^): the re- 
quirement that V is even preserves the parity rules used to go 
from dimerization to maps. 

The map uh factors as follows: 



riH : V^ 

I 

^3% 



(Tl) 



Uh 



We have checked this by (laboriously) producing the se- 
quence of plaquet flips carrying S to its refinement S' . But 
the vector fields ncont and n^cont associated to any dimeriza- 
tion and its refinement are (canonically) homotopic, as they 
always lie within 90° after scaling (xl/3) the refined lattice 
back to the original size. Thus Nh also factors 
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Nh: Vl 




> Z 



(T2) 



and by Theorem 3, Nh induces an isomorphism on ttq. So 
we have a diagram 



(T3) 



>Z2 




in which three triangles are known to commute: lower-left 
(Tl), upper-left (T2), and lower-right (a consequence of Nh 
being an isomorphism and fin an epimorphism). By abstract 
nonsense, the forth, upper-right triangle must also commute, 
which was to be proved. D 



D. Effective Field Theory for the 3D Dimer Model on the 
Cubic Lattice 

Let us now suppose that the dimers are governed by the 
Rokhsar-Kivelson Hamiltonian, 






i.c.\. 



(7) 



Then, it is known that fox V > t the system is in a stag- 
gered phase with no flippable plaquettes. At the RK point, 
t = V, the ground state can be found exactly: it is an equal 
amplitude superposition of all possible dimer configurations. 
For Vc < V < t, the system is in a Coulomb phase with 
power-law dimer-dimer correlations^. For V < Vc, the sys- 
tem spontaneously breaks translational symmetry and orders, 
depending on the value of V, in either a columnar, plaquette, 
or cube phase. 

We now reinterpret these results in light of our conclusion 
that a configuration of the 3D cubic lattice dimer model can 
be represented by a unit vector field n, at least at a course- 
grained level. We begin by writing the action for a continuum 
field theory description of the Rokhsar-Kivelson on the cubic 
lattice in the form S[n\. In the low-energy, long- wavelength 
limit, we can perform a gradient expansion for the action for 
this field: 



-I 



d^xdr 



^(a^n)2 + y(n) 



(8) 



not have 0(3) symmetry. Different possible forms for the po- 
tential V{n) will favor different possible ordered states. For 
instance, if we take 

V{n) = a{nl^n^^nj)^b{nl^n^^njf^cnlnlnl (9) 

then there will be regimes of the couplings a^b^c in which 
each of the three phases columnar, plaquette, and cube occurs. 
V{n) is clearly highly-relevant when ft orders since it deter- 
mines which ordered phase the system settles into. However, 
it may be irrelevant at the critical point to the Coulomb phase 
and it is certainly irrelevant in the Coulomb phase itself. 

It will prove to be very useful to use a different form for the 
action which is equivalent to ^ in the ordered phase. We will 
use the CP^ representation for the order parameter: 



(TxdT 



\{id^^a^)za 



\{zlZa - 1) 
-^V{Z^GZ)\ (10) 



In this equation, n = z'^gz. The coupling constant g con- 
trols the strength of fluctuations. For g large, there are large 
fluctuations in the direction ofn = z^az\forg small, fluctu- 
ations are small. The field A is a Lagrange multiplier which 
enforces the constraint z^Za = 1. The gauge field a^ has no 
kinetic term; it eliminates the phase degree of freedom of z, 
which does not enter into n. Hedgehogs in ft are magnetic 
monopoles of a^ : 

— / rnxe^:,x^ijknid^njd:,nk = ^ rnxe^^xf^iv (H) 



27r 



where /^^ = d^a^-dya^.a^ = ^{d^z^ z- z^duz),2indmx 
is the unit normal to the surface enclosing the monopole. 

In the ordered phase, (za) ^ 0. The gauge field, a^, is 
massive by the Anderson-Higgs mechanism. Consequently, 
monopoles/hedgehogs are confined. In the disordered phase, 
(za) = 0. The field Za can be integrated out, thereby gener- 
ating a Maxwell term for the gauge field a^. Thus, the disor- 
dered phase of this model is a Coulomb phase with a gapless 
photon. 

This is a slight refinement of the description of the U(l) 
gauge field description of the cubic lattice quantum dimer 
model^ which defines the magnetic field through: 



B^ir) = (-1) 



Qr 



r,r+i 



(12) 



where Q = (1/6, 1/6, 1/6), b is the lattice constant, and i = 
x^y^z. The Lagrangian for this field is, in the continuum limit, 

C = {dtSif - r(V X a)2 - k{V xV x a)^ (13) 

where B = V x a. At the Rokhsar-Kivelson point, t = V, 
r = 0. Otherwise, we can neglect the k term in the infrared 



Here, the potential V{n) is non-zero because the dimers do 



limit. The Lagrangian ( 13 ) is the low-energy limit of (10) in 
ao = gauge. The action ([TO]) also specifies the low-lying 
massive (in the Coulomb phase) degrees of freedom and gives 
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a picture for the gapped phases in terms of the condensation of 
these degrees of freedom. This description is also reminiscent 
(albeit in a different dimension) of the 2+1-D non-compact 
CP^ model of Motrunich and Vishwanath^, which describes 
the paramagnetic phase of an 0(3) a-model with hedgehogs 
suppressed. (See also Ref. |7l) 

A Hopfion configuration can be written in terms of the CP^ 
representation as a configuration of z which is topologically- 
equivalent to: 



2Arsin(9e^^ 



2 + A2 \2\t cos i9 + z(r^ - A^) 



(14) 



In the ordered phase, there is a gradient energy J ^x (pn)^ 
for such a configuration. Alternatively, we can describe this 

a^ which results when 



energy as the Meissner energy J d'^xa 



condenses, with a,. 



m 



r.t : 



z^dj^z). This energy 



scales as b if we rescale a Hopfion X ^ bX; therefore, a Hop- 
fion minimizes its energy by becoming as small as possible. 
The minimal-energy Hopfion will be only a few plaquettes in 
size - essentially a fluctuating version of the classical configu- 
ration depicted in Figure |3]- and its energy will be determined 
by lattice- scale dynamics. 

In the Coulomb phase, on the other hand, there is no stiff- 
ness for the field n. Instead, there is simply the Maxwell en- 
ergy 



U 



I 



d^x{E'^ ^B'^), 



which energy scales SiS U ^ U/b under a dilatation X ^ bX 
of the Hopfion ([14]). Therefore, an infinite-size Hopfion will 
have vanishing energy in the Coulomb phase of the 3D cubic 
lattice dimer model. 



III. FERMION DIMER MODEL 



In this section we define the fermion dimer model. We de- 



fine this model in Sections IIIB IIIC in a general fashion. 



such that the definition can be applied to any arbitrary lat- 
tice for which we can find a matrix M obeying the 7r-flux 
condition of Eq. ([T]). We then specialize to two dimensional 
and three dimensional examples in later subsections, show- 
ing how to realize deconfined Majorana excitations in three 
dimensions. 



A. Free Fermions on the Cubic Lattice with tt flux 

In Ref. 3 , we introduced a model of fermions hopping on 
a hypercubic lattice with tt flux through each plaquette. Spe- 
cializing to the case of three dimensions, there is a single Ma- 
jorana fermion operator 7r at each site r G Z^ of the lattice. 
We will have a unit cell with 8 sites, so this can be viewed 
as a model of spin- 1/2 electrons with a two-fold orbital de- 
generacy per unit cell (and charge non-conserving terms, such 
as would be present in a superconductor). The Hamiltonian 



takes the following form if the hopping has uniform magni- 
tude (with signs determined by the 7r-flux rule): 



H 



rGZ3 i=l 



I tr^r-|-x^7r7r-|-x^ 



(15) 



where t 



r,r+xi 



f f_nr-(x2+X3) f _ f f_nr-(x3) 

and tr,r+x3 = tz. Here, 7r is a real fermionic operator, 
7r = tI- This Hamiltonian is translationally-invariant since 
the flux, which is gauge-invariant, is the same through each 
plaquette. However, it is convenient to take an 8 -site unit cell 
on the cubic lattice so that the lattice constant is now 2. Then, 
we can group the 8 Majorana fermion operators in each unit 
cell, 7r+nx^, with r G (2Z)^, n = 0, 1, and i = 1, 2, 3, into a 
single 8-component spinor: 



Xa/37(r) = 7r+axi+/3x2+7X3 



(16) 



with a, /3, 7 = 0, 1. In the continuum limit, the Hamiltonian 
takes the form 



H = xocjidjX 



(17) 



where 



0-2 



I ^ax C 

/(g)/(g) 



(18) 



Here, we have specialized to the case tx = ty = tz, but taking 
anisotropic hoppings will simply make the velocity anistropic 
in the continuum limit, so long as the hopping strengths re- 



main non-zero. In Eq. 16 the first Pauli matrix acts on the 
index a in Eq. 16 the second Pauli matrix acts on /3, and the 
third on 7. 



we 


define 














>^U 
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1 

a/2 


(Xo/ii. 


+ Xi/ii.) 
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1 

V2 


K)m 


a (XOaz/ — 


X^lcxv 



then Eq. 17 can be re-written in the form: 

H = X^ (XjidjX^ + }? (XjidjX^ 



(19) 



(20) 



where ai = ciz^ (Tz, 0^2 = ctx^^-z, 0^3 = / (8) (Jx- From Eqs. 
16 19 and the reality of 7r, we see that A* = (A*)*. There- 



fore, (20) is the Hamiltonian of two 4-component Majorana 
spinors. In Eq. [TT] they have been been put together into the 



8-component field x defined in Eq. 16 

Equivalently, A^'^ form a single 4-component Dirac spinor: 

1 



V^ 



V2 



(A^+iA^) 



with Dirac Hamiltonian 



(21) 



(22) 
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The reason that we have introduced the 8-component field x 
defined in Eq. 16 and the Hamiltonian ^Tl) is that we will soon 
be considering mass terms (and, eventually, interactions with 
other fields) which are not invariant under the U(l) symmetry 



respected by Eq. 22 



Now suppose that the hopping strengths do not have uni- 
form magnitude. If they are staggered in the z-direction, 
so that tr,r+x3 = 1 + (— l)'^'^^ m while tr,r+x3 = 1 and 
^r,r+x3 = 1 are unchanged, then a mass gap opens. The 
Hamiltonian can be written in the form: 



H = xajidjX^rnxP3X 



(23) 



where ^s = il ^ I ^ cjy. If, instead, we had staggered the 
hopping strengths in the x-direction or ^-direction, we would 
have added a mass term with /3i = idy ® cFz ® cFz or P2 = 
il ^ (Ty^ (Jz, respectively. 

Suppose now that we have staggered the hopping strength 
in the z-direction by an amount m, as above. We decrease the 
hopping strength in the z-direction until we reach tz = m, at 
which point the hopping strengths in the z-direction alternate 
between 2m and 0. This does not close the gap. We can now 
reduce tx^ty, without closing the gap, until we reach tx = 
t^y = 0. Then, the system breaks up into an array of 2 -site 
systems, with a 2-level system i7r7r+x3 = ±1 for r • X3 G 2Z 
on each such 2-site systems (or for r • X3 G 2Z + 1 if we had 
taken tz = — m). We could similarly stagger the hopping 
strength in the x-direction and tune to the point tx = m^ty = 
tz = or stagger the hopping strength in the ^/-direction and 
tune to ty = m^tx = tz = 0. In these extreme limits, the links 
with non-zero hopping strength form a dimer configuration. 
The goal of this paper is to analyze what happens when these 
dimers become dynamical. 



where Uq » l;i ranges over the sites of the lattice; and 
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W, 



( Yl ^^j - 1) 



(26) 



<ij> 



Here, the sum over j ranges over sites which neighbor site 
i, and fiij is +1 if there is a dimer on the bond connecting 
site i to site j, and otherwise. Thus, Wi is diagonal in the 
dimer basis described above and equals if there is one dimer 
touching site i and is positive otherwise. 

The term Hff couples the fermions to the dimers. In the 
previous section, we defined a matrix M obeying the 7r-flux 
condition of Eq. ([T]). We use this matrix to define 



H 



ff 



= ^M, 



ij'^ij^ilj 1 



(27) 



where 7^ is the Majorana fermion operator on site i. 

Before defining the term H^^^, we pause to describe the 
physics of the Hamiltonian 



Hq = H. 



penalty 



H 



ff- 



(28) 



This Hamiltonian is diagonal in the dimer basis, so eigenstates 
of Hq can be chosen to be eigenstates of the dimer number op- 
erators fiij. For any given dimerization pattern (i.e., for any 
choice of whether each bond has a dimer or not), the eigen- 
states of Hq are 



\^d)^\^f{d)), 



(29) 



where ipd is a vector in the dimer Hilbert space with the given 
dimerization pattern d, and V^/(<i) is an eigenvector of the free 
fermi Hamiltonian 



Hff {d) = Y^ MijUij {d) ji-fj . 



(30) 



B. Non-Dynamical Dimers 



The Hilbert space of the system is the tensor product of a 
Majorana fermion Hilbert space and a dimer Hilbert space. 
The Majorana fermion Hilbert space has one Majorana mode 
per site of the lattice. The dimer Hilbert space has one two- 
state system on each bond. On each bond of the system, a 
dimer can be present or absent. Thus, the Hilbert space in- 
cludes states in which a given site has multiple dimers or zero 
dimers; however, we will add a penalty term to the Hamilto- 
nian to penalize these terms. States with distinct dimerization 
patterns are orthogonal. The Hamiltonian we study is 



H = H, 



penalty 



H 



ff 






(24) 



where i^penaity is a term penalizing any state which does not 
have one dimer per site. We choose 



H, 



penalty 



f/o^W^i, 



(25) 



where the numbers n^ (d) are equal to or 1 and are eigen- 
values of the corresponding operator fiij so that 



;#d) 



rii. 



■{dMd 



(31) 



For sufficiently large Uq, the ground states of Hq have 
dimerization patterns with no defects, so that each site has 
exactly one dimer touching it. To see this, note that for any 
dimerization pattern d, the smallest eigenvalue of Hff {d) is 
greater than or equal to — ^- nij{d), i.e. —1 times the num- 
ber of dimers in the given dimerization pattern. If the pattern 
has no defects, the smallest eigenvalue is exactly equal to —1 
times the number of defects, and the corresponding eigenvec- 
tor is: 

\i>°{d)} (32) 

which is defined by the condition 

Mijnij{d)jaj\i^°{d)) = -|V?(rf)) (33) 

for all i,j with riij = 1. In a dimerization pattern with de- 
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FIG. 14: Labelling of points in aplaquette. 



FIG. 13: Dimerization pattern with two defects with zero modes on 
each defect. 



fects at which more than one dimer touches a site, the small- 
est eigenvalue of Hff{d) is strictly greater than —1 times the 
number of dimers. Thus, while a dimerization pattern with 
defects can reduce the lowest eigenvalue of Hff{d) below 
—N/2, for sufficiently large Uq the penalty term for such pat- 
terns is larger than the reduction in fermionic energy and so 
for sufficiently large Uq, the ground states of Hq indeed have 
no defects. 

The Hamiltonian Hq has a highly degenerate ground state 
subspace, in one-to-one correspondence with the dimer con- 
figurations of the lattice. For large Uq, the low-energy excita- 



RK dimer modeP. It is equal to 



hL = vE{\\\) 



n 



(H) 



m 



' Swap]^3 



i.c.y (35) 



tions of Hq ( 28 ) above any of these ground states are gener- 
ated by acting with Majorana fermion operators. 



We include the "/" in the superscript of i^^^ to denote that 
this Hamiltonian couples the fermions and dimers, while in 
later discussions we use Hjix to denote the RK dimer model 
Hamiltonian. In the above equation, the sum ranges over all 
plaquettes of the lattice. The first term is the usual RK po- 
tential energy term. For V > 0, it penalizes configurations 
with parallel dimers on a given plaquette. The second term is 
very similar to the RK kinetic term, with one twist. It is the 
product of two different operators on the plaquette, one acting 
on the dimer Hilbert space and the other acting on the fermion 
Hilbert space. The first term, which acts on the dimer Hilbert 
space, does a plaquette move of the dimers. The term Swap^g 
acting on the fermion Hilbert space interchanges the fermions 
on sites 1 , 3 on opposite corners of the plaquette (we label the 
sites around the plaquette in clockwise order 1, 2, 3, 4 starting 



where m.i, = for all 1 < /c, / < p. Such a state has energy ^^ ^^^ ^^^ ^^f^ ^, ^^^^^ in Fig.fTg so that 
2p — ^ . • riij (d). Therefore, Hq has an energy gap equal to 2 



above the ground state subspace. 



7i -^ ±73, 
73 -^ T7i, 



(36) 



C. Dimer Dynamics 

The Hamiltonian Hq partly fulfills our goal, since 
monomers, which support Majorana zero modes, are decon- 
fined, as shown in Fig. [13] However, due to the enormous 



ground state degeneracy, this situation with deconfined de- 
fects is really a result of fine-tuning the Hamiltonian. Adding 
local interactions between dimers can lift the ground state de- 
generacy and produce interaction between defects. We now 
add such interactions between dimers by adding to the Hamil- 
tonian a term H^j^. Our goal is to construct H^j^ such that, 
for d = 3, the system has a unique ground state (or possibly 
an 0(1) or polynomial ground state degeneracy, depending on 
the topology of the lattice), with deconfined defects, and also 
such that the system is stable to adding additional weak inter- 
actions. The resulting Hamiltonian will have gapless gauge 
modes; in the discussion we consider the question of whether 
it is possible to construct a gapped system. 

The term HJ^j^ that we add is strongly reminiscent of the 



with the signs chosen so that if the configuration has no de- 
fects, and if the plaquette move changes the dimer configura- 
tion from dtod', then if the fermions are in the ground state of 
the Hamiltonian Hff{d) before the move, then they are in the 
ground state of Hff{d') after the move. To explicitly specify 
the signs, if the initial configuration has ni2 = 77,34 = 1, and 
ni4 = 7123 = 0, then the plaquette move turns this into a con- 
figuration with ni4 = 7^23 = 1 and 7112 = 7134 = 0, while the 
swap operator applies 



71 ^ (M34/Mi4)73 = -M34M1473, 

73 ^ (Mi2/M23)7i = -M12M2371. 



(37) 



±73 and 



Due to the tt-Aux rule ([T]), we indeed have 71 
73 -^ T7i- 

Note that i^^^ requires there to be 7r-flux in each plaque- 
tte in the matrix M because the interchange of two Majo- 
rana fermions requires one of them acquire a minus sign as 



in Eq. ( [36] ) in order to preserve fermion parity, respecting the 
superselection rules. 

Note also that the kinetic term in H^^^ does not change the 
total number of defects in the dimerization pattern. Further, in 
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the sector with no defects, the terms HJ^^^ and Hff commute. 
Thus, to determine the ground state of the Hamiltonian H = 
Hq^Hj^j^ in the no defect sector, we can directly use previous 
results on the quantum dimer model. If i/jq = ^^ A{d)\ilJd) 
is the ground state of the RK Hamiltonian with given t, /7o, 
where tpd is a state with given dimerization pattern and A{d) 
are complex amplitudes, then 



Y,Md)\^d)^W{d)). 



(38) 



is the ground state of H, where i/j^ (d) is the ground state of 
Hff{d). That is, the fermions "go along for the ride" as the 
dimer s move. 

Consider the RK point, t = F, at which the ground state of 
the fermion dimer model may be found exactly: 



\o)=N'Y,m^\^'f{d)). 



(39) 



Here, A/" is a normalization constant. This state has Hfdm 
eigenvalue —N/2. The ground state(s) is/are an equal ampli- 
tude superposition of all dimer configurations which can be 
obtained from a fixed one by acting with plaquette moves. In 
three dimensions, there will be different degenerate states cor- 
responding to different Hopf numbers, according to the dis- 
cussion of Section [ll] We view these as Hopfion excitations 
"above the ground state", albeit with zero energy. In addition, 
if the lattice has periodic boundary conditions, then there will 
be different winding number sectors which are not mixed by 
plaquette moves and, therefore, there will be different degen- 
erate ground states. 



Now consider the excited states 



|n,^2,. 



>=N^]^d')^\^%d)). 



d' 



(40) 



where the dimerizations d' all have m monomers at fixed 
sites n,i2, • • • ,^m- ^^K ^^^^ ^^^ move the monomers, 
and 1^1, ^2, • • • , ^m > is an eigenstate of Hfdm with energy 
—N/2 + m/7o. This energy does not depend on the distance 
between the monomers. Thus, the fermion dimer model has 
deconfined monomers at the RK point (as the ordinary dimer 
model does). Moreover, the states 7^^ • • • 7ip |n, ^2, • • • , ^m > 
also have identical energy —N/2 + mllQ. Thus, the decon- 
fined monomers support Majorana zero modes. 

We do not know exact fermionic excited states of the 
fermion dimer model, but we can consider the ansatz: 



i d 

= ^a,^\i^d)^\^'f{dy,i). 



H 



RK' 



(41) 



but it is an 



Such a state is not an eigenstate of Hq -\ 
eigenstate of Hq with eigenvalue —N/2 + 2. From the def- 
inition of Hj^j^, we see that the commutator [HJ^j^^ji] is a 
fermionic operator which acts only on the plaquettes neigh- 



boring site i. Thus, 



{ai\HFDM\ai) >-N/2^2-ct 



(42) 



where c is a constant which is roughly equal to the probabil- 
ity in the state |0) that an arbitrary plaquette has two parallel 
dimers. This implies that for t <C 1 (but V = t), there is still 
a gapped band of fermionic excitations even when the dimers 
resonate, albeit a gap which is reduced compared to Hq. Of 
course, it is possible that a completely different fermionic 
state has lower energy than the ansatz ( [42| ). However, this is 
unlikely since the equal-time fermion Green function is short- 
ranged: (0|7i7j|0) = unless i and j are on the same pla- 
quette. Furthermore, the existence of zero modes at defects 
implies that there is a fermionic gap in the bulk. Thus, it is 
likely that the fermions remain gapped at the RK point of the 
Hamiltonian Hq-\-HJ^j^. Indeed, this same argument says that 
the fermions remain gapped even away from the RK point, so 
long as Uq is sufficiently large and t is sufficiently small; this 
is important for the next two subsections where we discusses 
the dimer dynamics in liquid phases near the RK point. 

On the other hand, there is a gapless Hopfion at the RK 
point: it is the equal amplitude superposition over all dimer 
configurations with Nh = 1: 



\H)=Af 



d3NH{d) = 



m ® \i'°Ad)). 



(43) 



In the fermion dimer model, the uh = A/^7:/mod2 is simply the 
fermionic parity of the ground state of Hff{d). To see this, 
let us define the fermion parity as 



ieE 



(44) 



Here £^ is a sublattice of the original lattice which contains 
half of the sites; x is one of the basis vectors of the lattice; and 
£^, f are chosen so that if i e E then i -\- x ^ E is not. Then, 
in the ground state of Hff (d) defined by Eq. 
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(-1) 



Nf 



Pf(iV) 



(45) 



Therefore, at least at the RK point, the Hopfion is a gapless 
fermionic excitation. However, this gapless fermionic excita- 
tion does not affect the Majorana zero modes any more than 
the gapless bosonic mode which is present at the RK point 
does. In Section |Vj we will discuss whether this survives be- 
yond the RK point. 

Finally, we note that since the fermion parity is necessarily 
a conserved quantity, if the Z2 invariant Pf (A/") of the dimer 
configuration does not match the fermion parity, then the sys- 
tem cannot be in the ground state of Hff (d). Instead, an ener- 
getic price must be paid since the fermions must be in an ex- 
cited state for the given dimerization pattern. Thus, although 
it is possible to add terms to the ordinary dimer model which 
mix different Z2 sectors - namely, terms which move dimers 
around loops longer than a single plaquette - the Z2 invari- 
ant is protected by superselection rules in the fermion dimer 
model. 
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D. Two-dimensional Lattices 



E. 3D Cubic Lattice 



Ait = V, the ground state of Hrk is given, in any dimen- 
sion, by the equal amplitude superposition for all dimer con- 
figurations. As we have argued above, at this point, there are 
deconfined monomers which support Majorana zero modes 
while the fermions remain gapped. We now consider what 
happens when we move away from the RK point. In two di- 
mensions, for t > V or t < V, the dimers ordei^. In both 
cases, this leads to a confinement of the defects. Consider for 



example the configuration shown in Fig. 13 which displays 
columnar order as occurs for large negative V; there is a lin- 
ear confining potential between defects. Further, it has been 
argued that even the point t = V is very finely tuned^, and 
is unstable to adding other local perturbations to the Hamilto- 
nian. 



However, we can instead consider the model on a trian- 
gular lattice, where there is a liquid phase with deconfined 
monomers ^^ for a range of coupling near the so-called "RK 
point" t = V. In this case, we expect that the fermion dimer 
model will show deconfined Majorana zero modes near this 
point. In this regime, we expect that the system is described 
by a topological quantum field theory (TQFT). This theory 
should be a subtheory of Z2 <S) I sing, where the Z2 discrete 
gauge theory contains particles 1, e, m, em (identity, electric, 
magnetic, and product of electric and magnetic), and Ising 
contains the particles l,cr, ?/^. The present triangular lattice 
theory should be a subtheory of this, as it contains the mag- 
netic particle m, corresponding to Z2 vortices in the dimer 
configuration. It also contains the ifj particles (an excited state 
of the fermions). It also contains the particle ecr (a defect 
carries Z2 charge and a Majorana zero mode). So, the com- 
plete set of particles of the triangular lattice model should be 
1, m, ea^ emcr^ ijj^ mip. There is one subtlety in this set of par- 
ticles. In the ordinary Z2 discrete gauge theory of a dimer 
model, the particle e corresponds to a defect site. Given a 
defect site, there are two possible topological sectors: draw 
a line from the defect site to infinity, and consider the sec- 
tors with either an even or odd number of dimers crossing this 
line. The particle e corresponds to the equal superposition 
of those two sectors, and em corresponds to the superposi- 
tion of those two sectors with opposite signs. However, in the 
fermion dimer model, the sectors with even or odd dimer num- 
ber correspond to different Hopfion number (recall that for a 
planar lattice without holes, uh = 1 for all configurations, but 
the presence of the defect allows uh to assume either sign). 
While we cannot take a superposition of two configurations 
with different fermion parity, the presence of the a particle on 
the defect allows us to fix this problem, so that the particle ea 
corresponds to the equal amplitude superposition of the even 
and odd dimer sectors with the sign of the a on the defect 
changed in the odd sector. Similarly, ema is the opposite am- 
plitude superposition of those two sectors, again with the sign 
of the a changed in the odd sector. 



In the three dimensional cubic lattice, the situation is more 
interesting. It is believed^ that Hrk on a cubic lattice has 
a stable liquid phase for V slightly less than t. This phase 
has a gapless photon-like mode and is, therefore, called the 
Coulomb phase. Defects have a power-law interaction via this 
photon-like U{1) gauge field; this power-law interaction is at- 
tractive between defects on opposite sublattices, but decays as 
a power of distance between the defects, so that a single defect 
is deconfined: it can be taken arbitrarily far away from other 
defects while still having a bounded energy. 

In the fermion dimer model, Hfdm = Hq -\- HJ^j^, the 
situation is expected to be very similar. The gapped fermionic 
band is expected to remain gapped in the Coulomb phase, 
though the gap will be decreased by a term proportional to t, 
just as at the RK point. Similarly, monomers supporting Ma- 
jorana zero modes are expected to be deconfined. The main 
difference between the Coulomb phase and the RK point is 
that the interaction energy between two monomers at distance 
X is expected to be 



Vj 



RK 



duo d?q 

2^ (27r)2 



(const.) e*' 



5{x) 



(46) 



at the RK point (note that the (ao(g, 0) ao(— g, 0)) correla- 
tion function, which determines the interaction between static 
charges, is independent of q at zero frequency at the RK point, 
which is anisoptropic between space and time) and 
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(47) 



in the Coulomb phase. 

The Hopfion is a gapless excitation at the RK point. We 
do not have a direct calculation of the Hopfion energy in the 
Coulomb phase. However, we present effective field theories 
in Section |IV] which describe such a phase, and they indicate 
that the Hopfion has vanishing gap throughout the Coulomb 
phase. 

We now describe the statistics of the zero modes. The result 
of the study in Refs. 2 3 was that in the models considered 
there, there are two inequivalent ways to interchange a pair 
of monomers. Under the interchange, one of the Majorana 
modes changes sign, so that the interchange of defects on sites 
i, j leads to either 7^ -^ ^j , ^j -> -7^ or 7^ -^ -^j , ^j -> 7^. 
Further, it was found that it was possible to produce the trans- 
formation 7i -^ — 7i , 7j -^ — 7j , which does not exchange the 
defects but changes the sign of both of them, without moving 
the defects but by simply evolving the mass term under an ap- 
propriate trajectory. We will find very similar behavior in the 
fermion dimer model. 

To describe the statistics, we need to modify the Hamilto- 
nian H so that the defects become mobile. We can add terms 
to the Hamiltonian allowing processes such as those shown in 



Fig. 15 a,b), moving the defect and the dimer (we show only 
some of the possibilities in the figure, but there are others). 
The process described in the figure is a product of two 
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d) 



FIG. 15: a),b) Illustration of possible processes moving a defect and 
a dimer. The dashed arrow is used to indicate "before" and "after" 
configurations. In (a), the defect moves along a diagonal, while in 
(b) it moves along a straight line. c),d) How to compute the signs in 
the defect motion. The arrows on the dimers indicate the arrows as in 
Fig.[T] while the dashed arrow again indicates "before" and "after". 
Solid circles are used to indicate the defect site, while an open circle 
indicates that the defect site changes sign. 



terms. First, there is a term changing the dimer number from 
1 to on one bond and from to 1 on another bond. Second, 
there is a term interchanging the Majorana mode on the defect 
site and on one end of the dimer. Due to the superselection 
rules, this term must change the sign of one of the Majorana 

to 



modes. For example, using the site labelling of Fig. 14 



describe different sites, in Fig. [T5[a) the defect moves from 
site 1 to site 3 so we interchange 71 -^ ±73 , 73 -^ ^71 . 
We choose the sign in the interchange so that if fermions are 
in their ground state for the dimer configuration before the 
move, then they remain in their ground state after the move. In 



Fig. 15 c,d), we show two different cases leading to different 
signs; in one case the defect changes sign and in the other it 
does not. If the arrow convention is chosen following Fig. [T] 
then 15 c) and 15 ^d) correspond to different choices of the 
plaquette on which the defect is moving, as the arrow direc- 
tion depends on the plaquette. The sign also depends upon the 
dimer configuration within the given plaquette. Consider the 



process of Fig.[T5[c); if the defect had moved on the same pla- 
quette, but the initial configuration had instead a dimer con- 
necting the top two sites of the plaquette (rather than connect- 
ing the left two sites as shown), then the final configuration 
would instead have a dimer connecting the right two sites of 
the plaquette, and also the sign would be different: the defect 
would move without a sign change. Note that this dependence 
of sign upon dimer configuration only occurs when a defect 
moves diagonally within a plaquette, but not when it moves in 
a straight line as in Fig. (p3))(b), where the sign is fixed. 

We can also add terms to the Hamiltonian which favor the 
defect lying on certain sites. These are terms which reduce the 
penalty for having zero or more than one dimer touching a site. 
For example, we can modify Hpenaity which was previously 
defined to be Uq ^- Wi, so that it is instead ^^ Uo{i)Wi, al- 
lowing a site dependent penalty Uq. By reducing Uo{i) for 
certain sites i, we favor having the defect on those sites. 

By time-dependent control of these two terms, allowing de- 
fect motion and favoring certain sites for the defect, we can 
interchange two defects. That is, given defects at sites i, j, we 
can perform a sequence of transformations so that the defects 
end at sites j, i. Let us imagine that the terms in the Hamilto- 



nian are adjusted such that the path that the defect follows is 
almost completely deterministic: for a given defect, we pick a 
sequence of sites ^o, ^i, ^2, ••• We then require that the defect 
start at site i. We then add a term to the Hamiltonian allowing 
the defect to move between sites zq, H (which are assumed to 
lie a distance 2 apart so that this motion is possible) but not 
to any other sites, and gradually adjust the position-dependent 
Uo{i) so that the defect moves to sites n after a given pe- 
riod of time with probability close to 1. We then repeat the 
process moving from ii to 12, and so on. We say that this 
is "almost completely" deterministic because the defect may 
resonate back and forth between sites io and ii for some time, 
and so on, but otherwise the sequence of sites visited is com- 
pletely determined. We also assume that in this process we 
adjust the terms allowing the defect motion to only allow one 
of the two possible choices for a defect to move between two 
signs which are diagonally opposite on a plaquette (i.e., we 



allow either the process of Fig. 15 c) where the dimer moves 
from the left side of the plaquette to the bottom or the process 
where the dimer moves from the top side to the right side, but 
not both processes). 

We find that the sign a given defect acquires during its mo- 
tion is completely determined by the path it follows, i.e., the 
particular sequence of lattice sites it moves through, as well 
as the choice of dimer configurations in a plaquette when it 
moves between two sites diagonally opposite on a plaquette. 
As a result, there are, in principle, four distinct possibilities 
when we interchange defects: each defect can either acquire a 
minus sign or not. 

In two of the cases, when either both defects acquire a mi- 
nus sign or neither defect acquires a minus sign, a Hopfion is 
emitted. To understand this possibility, note that we can de- 
fine the Z2 invariant above even in a system with defects. We 
simply compute the Pfaffian of the submatrix of the matrix TV 
containing only the sites with no defects. If the initial config- 
uration has a given sign of the Pfaffian (the initial state may 
be a quantum superposition of different defect configurations 
with the given sign), then the superselection rules show that 
the sign of the Pfaffian changes after interchange. One way 
to prove this is to define the Z2 invariant for a system with 
defects in a slightly different way as follows: take the ma- 
trix N defined previously and then define a matrix N' which 
has the same matrix elements as N except if sites i, j are the 

two defect sites then N-, = ±i. Then, the matrix N^ has 

''J 

non-vanishing Pfaffian and its sign is an invariant under pla- 
quette flips. We choose some arbitrary choice of the sign of 
Nij for the two defect sites initially, calling one defect site 
i and the other site j and setting N-j = -\-i. Then, as the 
defects move, we define a sequence of different matrices N', 
one such matrix for each defect position, in the natural way: 
for each defect configuration with defects at sites /c, /, we pick 
the sign in the definition of Nj^i so that if the defect originally 
at site i moved to site k and changed sign a total of ni times 
along the process (recall that each time a defect moves, it may 
or may not change sign, depending on which sites it moves 
between) and the defect originally at site j moved to site / 
and changed sign a total of 712 times along the process, then 
Nl^i = +z(-l)^i+^2 Then, the sign of the Pfaffian of N' 
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does not change along this process. Thus, after interchang- 
ing the two defects, if neither one of them acquires a negative 
sign, then the sign of the Pfaffian of the submatrix of N' con- 
taining the signs without defects must change sign. That is, a 
Hopfion is emitted. In Section |VJ we will discuss its physical 
implications. 



IV. EFFECTIVE FIELD THEORIES OF THE COULOMB 
PHASE OF THE FERMION DIMER MODEL 

A. CP^ Model 

In this section, we give an effective field theory description 
of the fermion dimer model. We take the CP^ description of 
the ordinary cubic lattice dimer model, discussed in section 
as our starting point: 
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d^xdr 



Y K^^M + «/x)^a|^ + K^lZa - 1) 



■V{z^Bz)\ (48) 



We couple this Lagrangian to fermions on the cubic lattice 
with TT flux through each plaquette As discussed in Section 



III A the low-energy theory of such fermions is a single 8- 
component Majorana fermion: 



S = j d^xdr [x^drX + X^(^j^djx\ 



(49) 



where the a^s are defined in Eq. 



As discussed in Section III A 
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the effect of staggering the 
hopping strengths is to add mass terms to the Hamiltonian, as 
in Eq. ( [23] ). In the fermion dimer model, the dimer s effect an 
extreme form of staggering on the fermion hoppings. Thus, 
the natural coupling of the fermions to the dimers is: 



S = I d^x dr [x^drX + X^<^j^^jX + ^ ^^(^kZ xPkX 
1 



+ ^l(^^/i+^/i)^a| +^(^a^a-l) 



(50) 



where 
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= il ® I ®ay . 



(51) 



If the dimers order in, for instance, the columnar phase with 
the dimers aligned along the z-direction, then the fermion ac 
tion takes the form 



d^x dr [x^drX + X^(^j^djX + rn xp2.x] 



d xdr 



ip^drX + ip^ajidj^p -\- m ip^ Ps^p 



/ d^x dr [ilj^^id^i/j + m i/jiIj] 



(52) 



where ?/; is defined in Eqs. 19 and 21 , the ajS are defined after 

Eq. 20, ^s = I ^ ay; 70 = (3s, 7^ = ^sai, ip = ip^^o- In the 



/^s = / (8) cr^; 70 = Ps, 7i 

columnar phase, the action for the dimer degrees of freedom 
is most simply written in the form 



Sdh 



I 



d^xdr 



^g 



{d^fif + V{n) 



The interaction between the fermions and the dimers is: 



d^xdrrmlj^^ki^{nk - Sks) 



(53) 



(54) 



If the fermions are integrated out, the terms in 6'dimer are 
renormalized and a topological term is generated, as shown 
by Ran, Hosur, and VishwanatlP. This topological term takes 
the form: 



^topo — 

1 

19207r 



d^x dr e^^^^^Tr {-f^VdaVd(3Vd^Vd^Vdj,V) 

(55) 



where V takes values in SU(3)/SO(3) and is an extension of 
n to a 4 + 1 -dimensional manifold whose boundary is 3 + 1- 
dimensional spacetime; on this boundary, V takes values in 
5*^ C SU(3)/SO(3). There isn't a unique extension of n to V, 
but the resulting value for 5'topo (which must be an integral 
multiple of tt) is unique up to 27r. Therefore, e*'^*°p° is either 
+1 or —1, depending on whether two Hopfions are exchanged 
or not. 



In order to discuss the Coulomb phase in which n = z^az 
is disordered, it is useful to rotate the fermions to the local 
direction of the dimer order parameter, which is well-defined 
so long as hedgehogs don't proliferate. We rotate 



X^^X 
where IZ is defined by 

n^ {z'^akz) hn = h 

We can write IZ explicitly in the form: 

7^ = iii (/ (g) / (g) /) + iviMs - iu2M2 + iv2Mi 



(56) 



(57) 



(58) 



where Zj 


= Uj -\-iVj, and 




Ml = I <S)cry(^ax , 




M2 = (7y^(7z^ (Jx 




M3 = Gy^Gx® (Jz 



(59) 



Note that 7^ is a real matrix since the M^s are purely imag- 
inary; therefore, V^ = IZ^ . Then, in terms of the rotated 
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fermions, the action takes the form: 



S = d^xdr [x^drX 



+ X ajidjX^rnx f^sX 



X^n^drUx + x^^^o^jidjUx 

1 2 



1) 



(60) 



or, re-writing this in terms of the Dirac fermion ip, as in Eq. 
122 



S 






^\{id^^a^)Za\^ ^\{zlZa 



(61) 



where 1z = H^IZH and H 



1 

V2 



)/. 



We now consider the Coulomb phase, in which the dimers 
are disordered. In this phase Za is gapped, so we can inte- 
grate it out. This generates a Maxwell term for the gauge field 
and also a minimal coupling between the gauge field and the 
fermions: 



-I 



d^xdr 



-T , ^ 



V^7/i(^S/i + a^)^ + mV^V^ + -/^^ 



(62) 



?/; and a^ are minimally-coupled because the rotated ?/; is lin- 
ear in Za. Since Za is charged under the gauge field a^, the 
rotated ijj is as well. Therefore, the effective field theory for 
the liquid phase of the fermion dimer model is simply 3+1-D 
QED, but the fermions and the gauge fields may have different 
velocities (which will, of course, be different from the speed 
of light) and the fermion mass will be different from the bare 
electron mass. 

We note that we could have instead arrived at the action 
([62]) through a slave fermion construction. Suppose we define 
abosonic 'holon' Za and 'spinon' /q,/3^ through the definition: 



[(2:1 + zj) (/ (g) / (g) /) + {zi - zl)Ms 

■i{z2^zl)M2^{z2-zl)M^]f 



(63) 



This representation is redundant, as evinced by its invariance 
under 

/ -^ [cos 6> / / / + z sin 6> cr^ cr^; /] /. (64) 

This redundancy can be manifested, at low-energies, through 
the emergence of a U(l) gauge field a^ which couples to both 
Za and fo,j3^ . When Za is gapped, it can be integrated out, giv- 
ing Eq. 62 with ip replaced by the Dirac fermion constructed 



from / along the lines of Eqs. 19 and 21 



Note that this field theory describes the low-energy regime 
below the hedgehog (or, in QED language, monopole) gaps 
(assuming that the fermion gap m is smaller than the hedge- 



hog/monopole gap), where we were able to safely integrate 
out the field z. In order to discuss the physics of hedgehogs, 
we need to retain the field z since a hedgehog is configuration 
such as 



i0/2 



g-^|gi(/>/2 



(65) 



where r^O^cj) are polar coordinates. However, we will, instead, 
consider a different but related effective field theory in the 
next subsection, in which hedgehogs are explicitly retained, 
although the relation to the fermion dimer model is less di- 
rect. 

Hopfions, however, are non-singular configurations of the 
field z. Therefore, they are part of our theory, and are mani- 
fested as configurations of the gauge field a^, given by 



^/i = -{d^z^ z- z^dyz) 



(66) 



where z is a configuration which is topologically-equivalent 
to a Hopfion, such as: 



z = 



2Arsin6>e* 



+ A2 I 2Ar cos i9 + z(r2 - A^) 



(67) 



As noted in Section II D such a configuration has vanishing 
energy as A -^ 00. However, this gapless excitation is actu- 



ally a fermion, as a result of the topological term ( [55] ) which 
results when the gapped fermions are integrated out. There- 
fore, the action ( [62] ) is misleading; the system actually has a 
gapless fermion - the Hopfion. In Section [Vj we will discuss 
the effect of a gapless Hopfion. First, however, we will dis- 
cuss an effective theory in which hedgehogs and Hopfions are 
retained explicitly. 



B. 4 + 1 -Dimensional Model 

In this subsection we will propose an alternative effec- 
tive theory description of the Coulomb phase, which has a 
completely different starting point but leads to results con- 
sistent with the one proposed in the last subsection. This 
approach is motivated by the t'Hooft-Polyakov approach to 
magnetic monopoles^^ ^^. We start from the action of Majo- 
rana fermions with a mass term: 

S = j d^xdr [x^drX + X^(^j^djX + mnkXPkX 



1 
'2^ 



d^Ukd^Uk 



(68) 



where Uk is the order parameter ft written in components. It 
is natural to ask whether one can modify the model by in- 
troducing an SU(2) gauge field, which is broken to /7(1) by 
the mass terms. If this can be done, the linearly divergent 
monopole energy will be screened and the interaction between 
monopoles assumes a Coulomb form, similar to the t'Hooft- 
Polyakov monopole. Such an approach turns out to be impos- 
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sible due to an anomaly, as noted in Ref . T, but the problem 
can be cured by introducing a (4 + 1) — (i theory, as we now 
discuss. 

To introduce an SU(2) gauge field it is convenient to rewrite 
the action in terms of Weyl fermions: 

S = j d^xdr [c^ {dr + (J'i-idi)) c 

+ (mnkC^ {dy (g) TyTk) c^ + h.c.j + ^d^rikd^nk 

(69) 

in which c is a two-component spinor in space-time, and also 
carries an isospin 1/2 representation of an internal SU(2) 
symmetry. In components, one can write c^rs with a = 1, 2 
spin indices and 5 = 1,2 isospin indices, a* and r^ stands for 
Pauli matrices in spin and isospin indices, respectively. Phys- 
ically, one can consider this action as a spin singlet supercon- 
ductor formed by Weyl fermions. Due to fermionic statistics 
a spin singlet pair must be triplet in isospin, which is why n^ 
carries a vector representation of the isospin symmetry. From 
this expression, one can see that the kinetic energy terms are 
SU(2) invariant, and the mass term breaks SU(2) symmetry 
to U(l). In this theory the monopole has linearly divergent 
energy. To avoid such a divergent energy one can gauge the 
theory by introducing an SU(2) gauge field a^ : 



S' 



-I 



d^xdr \c^ (dr — iao + a\—idi — ai)) c 
- (mnkC^ {ay (g) TyTk) c^ + h.c.j 
-^ D^n^D^nk] + ^Tr [UJn (70) 



In such an action, the order parameter Uk carries gauge charge 
of the gauge field a^. When the order parameter rik con- 
denses, it induces a mass term for both the fermions and, 
by the Higgs mechanism, for the gauge field as well. The 
monopoles become U(l) magnetic monopoles of the resid- 
ual U(l) gauge field, so the interaction between them is a 
Coulomb interaction. 

However, the above reasoning fails when we take into ac- 
count the axial anomaly of the Weyl fermion. If we introduce 
an ultraviolet regularization of the theory above, the regular- 
ization has to break the SU(2) symmetry. More explicitly, if 
we consider a lattice regularization of the theory, according 
to the the Nielson-Ninomiya theorem^ ^ there must be dou- 
bling partners c of the Weyl fermions c with opposite chiral- 
ity. The UV regularization must break chiral symmetry so that 
the Weyl fermions and their doubling partners are coupled and 
there is only one SU(2) symmetry transforming both of them. 
To make the doubling partners gapped, there must be a mass 
term hk for c which breaks SU(2) symmetry. Thus the ef- 
fective action of the order parameter Uk in general contains a 
coupling term fikfik other than the ordinary kinetic energy and 
potential energy terms. With n a constant field determined by 
the cut-off, the monopole of ns will reobtain a linear energy 
cost, so that we conclude that in a properly regularized version 



of theory (70), the monopoles remain confined. 

Interestingly, a resolution to this problem can be found by 
introducing a (4+1) -dimensional theory. It is well-known that 
(3 + 1) — d Weyl fermions can be regularized as "domain wall 
fermions" living on the boundary of a (4 + 1) — d lattice Dirac 
theorylS. More generically, it was demonstrated in Ref.^^ that 
any (4 + l)-d gapped fermion theory with a nontrivial second 
Chern number in the geometrical gauge field defined in mo- 
mentum space has such domain wall fermions on its bound- 
ary, and the number and chirality of the Weyl fermions on the 
boundary is determined by the second Chern number. Such 
a (4 + l)-d state is a topological insulator which is the par- 
ent state for (3 + 1) and (2 + 1) dimensional time-reversal 
invariant topological insulators. ^^ 

Here we consider such a lattice regularization by using the 
lattice Dirac model. In real space 



54+1 = dr 
Jo 



i,Q;=l,2,3,4 



i i 

tir^ + ^r^ 

' — Ci^a + h.C. 



Cj 



(71) 



with r0'i'-'4 



Hermitian Dirac T matrices F^^ 



5^^. In the continuum limit this ac- 
tion describes a massive Dirac fermion with mass M. For 
the mass range — 2 < M/B < the Dirac model has Chern 
number C2 = 1 in the Brillouin zone, and on the boundary 
there is one Weyl fermion. To obtain the Weyl fermion SU(2) 



doublet in Eq. ( [7Q| ) we consider an SU(2) doublet of q in 
(4 + l)-d, coupled with an SU(2) lattice gauge field and a 
Higgs field: 






4+1 



/3 



dr 



Y^ 4 {dr - iao) Q + (M + 4B) ^ cJf^q 



2^ I q 7. e* ^"Q+a + h.C. 



;,a=l,2,3,4 



- ^ (mnkcl {T ryrk) c\ + h.c. 



(72) 



with SmI^h] ^^^ Maxwell term for a^ and 6'sigma[n] the lat- 
tice version of the sigma-model action ^D^UkD^Uk. T is 
the time-reversal matrix similar to ay for the two-component 
fermion. In the same way as in the (3 + l)-d theory, the order 
parameter n^ is a Higgs field carrying a triplet representation 
ofSU(2). 

Consider this theory defined on a finite thickness slab 
R^ X I with finite thickness along the 4-th dimension and in- 
finite along the other three dimensions. If n^ = 0, a doublet 
of massless Weyl fermions live on each of the two b ound- 
ary R^ with opposite chirality, as is illustrated in Fig. IV B 



(a). We can consider the Weyl fermions on the two bound- 
aries as the c fermions in Eq. ([7Q|) and their doubling partners 



20 



(b) 



^^ 




FIG. 16: (a) Illustration of the (4 + l)-d theory, which consists of 
lattice Dirac fermions coupled with a SU(2) gauge field described 
by the action ^. For -2 < M/B < the bulk is gapped 
and topologically nontrivial, with one doublet of Weyl fermions on 
each surface, (b) When the Higgs field ft is in the condensed phase, 
the surface Weyl fermions become gapped Majorana fermions, and 
the monopoles on the 3d surface are end-points of one-dimensional 
monopole lines in the 4d bulk. 



c. The mass term Uk breaks SU(2) symmetry and transforms 
the Weyl fermions into massive Majorana fermions. Since we 
consider the top surface as the physical system and the bot- 
tom surface as the doubling partners, we always consider the 
n configurations with fixed boundary condition n = no at 
the bottom surface, so that the UV cut-off is fixed and has no 
topological configuration. 

The advantage of the (4 + l)-d theory is that it describes 
the anomaly of the boundary theory explicitly. In (3 + l)-d, 
the point-like monopole of the n field can be defined since 
7T2{S'^) = Z. In (4 + l)-d the same homotopy group leads 
to one-dimensional topological defects. Thus the monopole 
on the boundary 3d space is always the end of Id "monopole 
lines" in the bulk, as shown in Fig. 



IV B (b). This picture 



immediately leads to the conclusion we obtained earlier that 
the monopoles are still linearly confined after the theory is 
gauged. Indeed, the linear confinement is due to the existence 
of the UV cut-off on the other boundary, since the boundary 
condition on the bottom surface does not allow the monopole 
lines to penetrate in the other surface. (Otherwise the lowest 
energy monopole lines should be perpendicular penetrating 
both surfaces when the monopoles are far away.) 

The existence of non-Abelian statistics can also be under- 
stood explicitly in the (4 + 1) dimensional theory. As is dis- 
cussed in Ref , 3 and earlier in this draft, the essential fact en- 
abling the possibility of non-Abelian statistics in (3 + l)-d is 
the existence of two topologically distinct exchange process 
of two monopoles, which are different by a "braidless opera- 
tion", i.e. spinning one monopole by 27r while keep the other 
monopole fixed. If we create a monopole-anti-monopole pair 
and do the braidless operation, and then annihilate the pair, 
the configuration left is not the vacuum but the Hopf map. In 
the (4 + l)-d theory, after such a process the ft configuration 
on the top surface is a nontrivial map with winding number 
1 in 7r3(6'^) = Z. In the bulk, such a process will create a 
point-like defect in the 4d space, just as 7r2(6'^) = Z leads 
to stable point-like defects in 3d space. Thus the "Hopfion" 
we discussed earlier is a point-like defect in the (4 + l)-d 
theory, which can be some distance away from our physi- 
cal 3d space on the boundary. This argument is illustrated 
in Fig. [17] In relation to the chiraly anomaly, it can be 
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FIG. 17: Illustration of the creation of a Hopfion in the bulk by creat- 
ing a monopole-anti-monopole pair, rotate the monopole by 27r and 
annihilate the pair. The Hopfion is a point-like defect of the order 
parameter n in the bulk 4d space. 



tion of the SU(2) gauge field with nontrivial instanton number 
3^ J d^xdre^^^^ ill f^j^f err] = 1, which thus pumps one 
fermion from the boundary into the bulk. After the monopole- 
anti-monopole pair is annihilated, the fermion pumped by this 
process stays with the Hopfion in the bulk. 

Now we address the question of whether it is possible to 
have a phase with deconfined monopoles in this theory. In the 
(4 + l)-d theory, the confinement of monopoles is due to the 
monopole line connecting two monopoles. This is the same as 
confinement between U(l) magnetic monopoles in a (3 + l)-d 
superconductor. To make the monopoles deconfined we need 
to suppress the Higgs phase. In last subsection, we quantum 
disordered the order parameter field z for this purpose. In 
the (4 + l)-d theory, the deconfined phase can be introduced 
in a different, somewhat more explicit way by making use of 
the extra dimension. To this end, consider a modified ver- 
sion of the (4 + l)-d theory as shown in Fig. 18 with the 
(4 + l)-d slab geometry the same as that discussed above but 
the order parameter n = in a region < X4 < a includ- 



ing the bottom surface and its neighborhood. The region with 
n = is a normal topological insulator phase (denoted by N) 
preserving SU(2) symmetry, and the rest of the system (de- 
noted by S) which is a "color superconductor" of the SU(2) 
isospin. The bottom surface state has a width of ^ ~ l/l^l 
in the x^ direction. If a :^ 1/|^| the bottom surface can be 
considered as decoupled from the symmetry breaking order 
parameter field ft and thus remains massless Weyl fermions. 
In such a geometry, a monopole on the top surface can be 
attached to a perpendicular monopole line along the x^ direc- 
tion which terminates at the domain between the N and S re- 



shown that the braidless operation corresponds to a configura- 



gions, as shown in Fig. 18 The interaction between two such 
monopole lines perpendicular to the (3 + l)-d boundary is the 
(3 + l)-d Coulomb interaction, just as for ordinary t'Hooft- 
Polyakov monopoles. In the N region, the SU(2) gauge in- 
variance is recovered and the SU(2) magnetic flux can spread 
in space. It is important to note that a Majorana zero mode 
is trapped at the end point of the monopole line at the top 
surface, but no Majorana zero mode exists at the bottom end 
point of the monopole line, since the fermion degree of free- 
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dom is completely gapped at the N/S boundary. Now we con- 
sider a rotation of the monopole on the top surface by 27r. If 
the bottom end point of the monopole line is not rotated, a 
Hopfion will be created in the S region which is a point-like 
topological defect costing some energy. If the N region is in 
the Coulomb phase of the SU(2) gauge field, the Hopfion will 
prefer to move down to the N region where it will smear in the 
whole space to save energy. As is discussed above, the Hop- 
fion corresponds to a soliton configuration of the SU(2) gauge 
field with nontrivial second Chern number, which thus carries 
one fermion due to the nontrivial response of the (4 + l)-d 
topological insulator. Because the bottom surface is now gap- 
less, if we rotate the monopole slowly the Hopfion will finally 
propagate to the bottom surface and disappear, so that the en- 
ergy cost from the gauge field strength and order parameter 
can be completely eliminated, and one fermion is pumped 
from the top surface to the gapless bottom surface. Thus, 
in this configuration we have realized deconfined monopoles 
with a single Majorana zero mode and Coulomb interaction. 
One can see that gapless Weyl fermions on the bottom surface 
are necessary for realizing this phase while still preserving the 
non-Abelian statistics of monopoles, without which it would 
be impossible to pump a fermion away from the top surface. 
For a slab with finite thickness L there is a nonzero coupling 
between the top and bottom surface states. Consequently, in 
principle the gapless bottom surface states will mediate a cou- 
pling between two Majorana zero modes on the top surface 
which decays as a power law rather than exponentially in dis- 
tance. However, such a coupling is suppressed by an expo- 
nential factor of e~^^ which can be made arbitrarily small 
for L > ^ = 1/M. Compared with the (3 + l)-d theory 
discussed earlier, we see that the Weyl fermions on the bot- 
tom surface in the (4 + l)-d theory correspond to Hopfions 
in the (3 + l)-d field theory, or the lattice fermionic dimer 
model. The (4 + l)-d and (3 + l)-d theories lead to con- 
sistent results for a deconfined phase with gapless fermionic 
Hopfions. Moreover, the (4 + l)-d theory provides an explicit 
way to describe the fermionic Hopfions and a controlled way 
to suppress the coupling between Majorana zero modes in the 
monopoles and the gapless Hopfions. 

It should be noticed that the discussion above can be gener- 
alized to more generic systems with A^ Weyl fermions coupled 
to an U(N) gauge field. The (3 + 1) dimensional action (70) 
can be generalized to 



S' 



d^xdr [c^ [dr — icio + cr^{—idi — a^)) c 
+ ^ Tr {D^n^D^n)] + ^Tr [f^.r] (73) 



with c a N-dimensional spinor carrying fundamental repre- 
sentation of U(N), Uah Sin N X N mass matrix satisfying n = 
n^, n^n = 1. This mass term gives the Majorana fermion a 
mass m and breaks SU(N) to the subgroup of SO(N). The or- 
der parameter manifold is thus SU{N)/SO{N). This model 
can be considered as a generic model for Majorana fermions, 
while the SU(2) model discussed above corresponds to the 
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FIG. 18: Illustration of the configuration with the order parameter 
in the Higgs phase in the S ("superconducting") region and vanishes 
in the N ("normal") region. Each perpendicular monopole line traps 
a single Majorana zero mode confined at the top endpoint, and two 
monopole lines are coupled by Coulomb interaction. The bottom 
surface state is a SU(2) doublet of gapless Weyl fermion. 



special case of A" = 2. In SU(2) case we have the Hopf map 
coming from nzi^S^) = Z. In the case of SU(N) (A^ > 2) 
we have 7r3{SU{N)/SO{N)) = Z2 so that only a Z2 non- 
trivial generalization of the Hopf map can be defined. We also 
have 7r2{SU{N)/SO{N)) = Z2 so that monopoles still ex- 
ist but there is no distinct monopole and anti-monopole. For 
such a model, the discussion above all holds since we can 
still define the (4 + 1) -dimensional model, which still have 
monopole lines and point- like Hopfons. By the same con- 
figuration as is discussed in Fig. 



18 a phase with weakly- 



interacting monopoles can be found. 



V. MONOPOLE STATISTICS 

We now discuss what happens when monopoles are ex- 
changed and twisted. (We refer to any process which changes 
the sign of a single Majorana zero mode operator and emits a 
Hopfion as a "twist" since it produces the same effect on the 
defect as the 27r rotation in Ref. 3 ) As discussed in the context 
of the dimer model, when two defects are interchanged, we 
have four possible choices: either, neither, or both Majorana 
zero modes could change sign. In the ordered phase (or any 
other phase in which a Hopfion costs non-zero energy), one of 
the zero modes must change sign (if the process is done adi- 
abatically); which zero mode changes sign is determined by 
the path. If a single monopole is twisted, then a Hopfion must 



be emitted, as we discussed in Section |IirE| (twisting a single 
monopole changes the sign of its Majorana zero mode so a 
Hopfion is created to preserve fermion parity). In the ordered 
phase, the emitted Hopfion must be virtual if the twisting is 
done slowly enough. Therefore, the virtual Hopfion must be 
absorbed by another monopole, which also twists, as a result. 
This is a process that changes the sign of both defects. 

In controlled processes, we hold the monopoles fixed and 
determine which monopoles twist. However, suppose that we 
allow quantum fluctuations of the defects in the fermion dimer 
model. (The Coulomb phase will generically have such fluc- 
tuations.) That is, suppose that we move defects with a posi- 
tion dependent Uo{i), but we have non-vanishing terms in the 
Hamiltonian allowing the defect to move locally. Then, an in- 
teresting process is possible in which the defect moves around 
a short loop, changing its sign and emitting a Hopfion. This 
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Hopfion can then be absorbed by the other defect, changing 
its sign. That is, there is an ampHtude for a process that spon- 
taneously changes the sign of both defects. When a Hopfion 
is transmitted from one defect to the other, it induces an effec- 
tive interaction term ^7^7^ in the Hamiltonian, where 7^,7^ 
are the operators corresponding to the zero modes on defect i 
and j and where the coefficient of this effective term depends 
upon the spacing between the defects. We expect that if the 
defects are far separated, then this process will be suppressed 
exponentially in the spacing between the defects. 

Now consider the Coulomb phase. The first question is 
whether the Majorana zero modes survive at all in the pres- 
ence of gapless Hopfions. At the RK point of the fermion 
dimer model, the zero modes are completely unaffected by the 



Hopfion, as discussed in Section [IirC| However, this is a spe- 
cial feature of the fermion dimer model. At a generic point in 
the Coulomb phase, we expect a coupling between the Majo- 
rana zero modes and Hopfions. However, this coupling need 
not destroy the Majorana zero modes. Consider the follow- 
ing toy model. Suppose we have a system of non-interacting 
Majorana fermions, arranged on a d-dimensional lattice with 
TT-flux hoppings. Add an additional defect site which couples 
only to a single site, which we call site 0, in the lattice, so 
that the full Hamiltonian is H = Hq -\- itjdefectjo, where Hq 
is the TT-flux Hamiltonian, t is the tunneling strength between 
the defect site and site 0, and ^defect and 70 are Majorana op- 
erators on the respective sites. This system has an exact zero 
energy state (exact because we added only one defect). The 
amplitude of this state on site x is proportional to Gqx^ where 
the Green's function G is G = H^^ . For d = 2, Gqx is 
proportional to l/\x\, so that the norm square is propotional 
to l/|xp and the integrated norm squared diverges logarith- 
mically; this implies that in d = 2 the state is delocalized. 
However, ford = 3, Gqx oc l/|xp, so the integrated norm 
squared converges and we have an exactly zero energy state 
localized near the defect. So, for d > 2, there is a local- 
ized state with an infinite lifetime, despite the coupling to a 
gapless continuum. Based on this toy model calculation, it 
seems likely that similarly the defect mode can have an infi- 
nite lifetime, despite coupling to gapless Hopfions, for d = 3, 
although some caution is required since this answer is clearly 
sensitive to the density of states of the gapless continuum near 
zero energy, and the answer could change for different disper- 
sion relations. 

We can analyze the effect of the coupling of a zero mode to 
a continuum more generally by considering an effective field 
theory of an isolated monopole at the origin with its associated 
zero mode. Let us suppose that the low-energy effective action 
for the Hopfion can be written as a relativistic Dirac fermion 
as in the toy model above. Then the effective action for the 
Hopfion and the Majorana zero mode will take the form: 

+ V / dTj{uilj{0) + c.c.) (74) 
where ^|J is the Hopfion annihilation operator, 7 is the Ma- 



jorana zero mode, v is the magnitude of the matrix element 
between the zero mode and the Hopfion, and u is 3. constant 
unit spinor which is the matrix element between the zero mode 
and the Hopfion divided by v. Then, by power counting, v has 
scaling dimension 1 — f . In three dimensions, the coupling 
between the zero mode and the bulk fermion is irrelevant and, 
therefore, unimportant in the infrared. The behavior seen at 
the RK point of the fermion dimer model is the fixed point 
behavior of the phase. This can also be seen by an explicit 
calculation of the zero mode Green function at lowest order in 



v: 



ili^^nhi-i^n)) = 



ICr, 



(75) 



where Z~^ = 1 + O(v^). There is still a pole at zero en- 
ergy. The zero mode is spread out over a few lattice sites due 
to its interaction with the bulk, so that the probability for it 
to remain at the origin is reduced by a factor of Z. How- 
ever, it has not disappeared into the bulk continuum, as would 
have happened if there had been a constant Hopfion density- 
of- states at zero energy, in which case the zero mode Green 
function would have had a form ^ (ze^ + iTsgn{en))~^ for 
some constant F. In the field theories of the previous section, 
the Hopfion has a relativistic spectrum and hence has vanish- 
ing density of states at zero energy; it is clearly not a Dirac 
particle, since it has additional degrees of freedom (a Hopfion 
can stretch and distort in various ways, as well as changing its 
overall size), but this actually helps us since the low energy 
Hopfions are necessarily large, so they will have smaller ma- 
trix elements to couple to a point-like defect. Thus, we expect 
that the zero mode survives in the full theory. 



In Eq. 74 we have assumed that Hopfion number is con- 
served, as it is in the fermion dimer model. However, if dimer 
moves involving more than one plaquete are included in the 
Hamiltonian, then the Hopfion number will not be conserved; 
it will only be conserved modulo 2, as is required by fermion 
parity. Thus, the generic situation is, in fact, that the Hopfion 
is a Majorana fermion. 

In the Coulomb phase, the amplitude for a virtual Hopfion 
to be exchanged between two monopoles will fall off poly- 
nomially with the distance between the defects. Therefore, 
the splitting between the two states of a pair of monopoles 
will fall off as an inverse power of the distance between them, 
rather than exponentially, as we would have in a true topolog- 
ical phase. For d = 3, using the non-interacting toy model 
above, the splitting is inversely proportional to the square of 
the distance. 

The gaplessness of the Hopfion has another important con- 
sequence. A single monopole can twist and emit a zero en- 
ergy Hopfion. If T] is the creation operator for a zero-energy 
Hopfion (which has infinite size since the Hopfion energy is 
inversely proportional to its linear size) and 7^ is the Majorana 
fermion operator for the i^^ quasiparticle, then twisting the i^^ 
quasiparticle implements the transformation 



Ri = r]li 



(76) 



Meanwhile, exchanging quasiparticles i and j without twist- 
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ing either one (which would be necessary to heal the order pa- 
rameter in the confining phase of the hedgehogs) implements 
the transformation 



Thus, the combination of an exchange and a twist is 



T,, 



-^i^ij 



1 



(l + 7z7,) 



a^lilj/^ 



(11) 



(78) 



which is the representative of a generator of the ribbon permu- 
tation group. Thus, the set of low-energy processes which are 
possible in the Coulomb phase is larger than simply the ribbon 
permutation group. Hopfion emission processes are possible; 
the ribbon permutation group is the subset of low-energy pro- 
cesses in which real Hopfion creation or annihilation does not 
occur. 



VI. DISCUSSION 

When a quasiparticle in a 21) topological phase is rotated 
by 27r, or twisted, the state of the system changes by a phase 
which depends only on the particle type. If there are several 
particles of the same species, then the phase acquired by the 
system doesn't depend on which particle is rotated. This fol- 
lows from locality: an operation performed on a single parti- 
cle cannot affect distant particles and must, therefore, be just 
a phase. However, in a system of 3D non-Abelian anyons 
- or, to be more precise, particles which obey projective rib- 
bon permutation statistics - a twist has a highly non-trivial ef- 
fect. Rotating a particle anti-commutes with rotating a differ- 
ent one. In the model introduced by Teo and Kane^ and in the 
general framework introduced in Ref . 3 , there are confining 
long-ranged interactions V{r) ex r between the hedgehogs 
which are 3D non-Abelian anyons. Thus, it is plausible that 
an operation performed on a single particle can affect distant 
particles. Indeed, the order parameter configuration which in- 
terpolates between hedgehogs can only be healed after an ex- 
change if one of the hedgehogs is twisted. Consequently, ex- 
changes belong to the ribbon permutation group^, rather than 
the ordinary permutation group. However, this alone does not 
explain why twists anti-commute. It must, furthermore, be 
the case - as it is in the model introduced by Teo and Kane^l- 
that the ribbon permutation group is represented projectively, 
rather than linearly. The multiplication rule of the ribbon per- 
mutation group is only represented up to a sign (but hardly an 
innocuous sign!) on the Hilbert space of the theory. This is 
presumably related to the fact that the overall phase associated 
with any exchange is not a topological quantity - since, as a 
result of the long-ranged force between hedgehogs, it depends 
on non-universal details of the exchange. Therefore, it need 
not linearly represent the ribbon permutation group. Thus, the 
fact that hedgehogs are not quite point-like quasiparticles is 
doubly important: (1) since there is an order parameter inter- 
polating between them, their exchanges belong to the ribbon 
permutation group; and (2) since there is a long-ranged force 



between hedgehogs, the ribbon permutation group is projec- 
tively represented. 

Thus, it is a little difficult to imagine how 3D non-Abelian 
anyons could be the weakly-coupled point-like quasiparticles 
of a system. Nevertheless, in this paper, we have given two 
complementary field theory descriptions of a Coulomb phase 
of 3D non-Abelian anyons. In this phase, twists of different 
quasiparticles do not commute because twisting a quasipar- 
ticle causes a Hopfion, a gapless fermionic excitation of the 
system, to be emitted. In fact, two twists anti-commute, as 



may be seen from Eq. 76 In this way, quasiparticles com- 
municate their twists to distant particles through the emission 
and absorption of gapless Hopfions. The very gaplessness of 
the Hopfions which facilitates the non-trivial statistics of the 
monopoles also reduces the protection of the quantum infor- 
mation which is contained in their zero modes. The splitting 
between zero modes falls off only as a power law, rather than 
an exponential. It does so as a result of Hopfion phase space 
restrictions, but this is no worse than a Fermi liquid, which is 
also stable due to phase space restrictions, rather than a gap. 

In a recent preprint, McGreevy and Swingle ^^ argue that a 
magnetic monopole cannot have a single Majorana fermion 
zero mode in a spontaneously-broken SU(2) gauge theory. 
Our phase does not conflict with their conclusions because 
they demand that the fermions be fully-gapped in the bulk 
while our model necessarily has a gapless fermionic Hopfion. 

We have also constructed an explicit microscopic model 
which has both confining and Coulomb phases of hedgehogs 
supporting Majorana zero modes, the fermion dimer model. 
The field theories constructed in this paper appear to embody 
the physics of the fermion dimer model. However, it is possi- 
ble, in principle, that the Coulomb phase of these field theo- 
ries is different than the Coulomb phase of the fermion dimer 
model. To show that they are the same, we would need to com- 
pute the energy of the Hopfion in the fermion dimer model. 
At the RK point, the Hopfion has vanishing energy gap. We 
would like to know if the Hopfion gap remains zero in the 
Coulomb phase. This is already an interesting question in the 
ordinary dimer model (i.e. without fermions), and a quantum 
Monte Carlo simulation could determine the energy gap of the 
Hopfion and study its behavior near the critical point. 

Can this model be used to perform topologically-protected 
quantum computation? Since the Hopfion is gapless, it re- 
quires a long time to perform an operation such as twisting 
both quasiparticles while staying in the adiabatic regime in 
order to return to a state without a Hopfion at the end of the 
process. That is, even though neither particle moves a large 
distance, the operation requires a long time; this should be no 
surprise given that the Hopfion must be emitted at one particle 
and absorbed at the other a long distance away. This long time 
poses a problem since we have also argued that if defects have 
a quantum dynamics, then they can absorb and emit virtual 
Hopfions leading to a weak level splitting; if the time required 
to perform a twist operation is sufficiently long compared to 
the inverse level splitting then the desired operation will not 
be impliemented with high fidelity. One way around this is to 
imagine that we "pin down" the quasiparticles and stop them 
from twisting; that is, we imagine that we perfectly control 
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the path along which each defect moves. 

The gapless Hopfion raises the natural question of whether 
we can construct a gapped Hamiltonian with deconfined non- 
Abelian anyons in three dimensions - i.e. a phase in which 
non-Abelian anyons have exponentially-decaying, rather than 
confining or even Coulomb interactions. At a field-theoretic 
level, one possible route is the condensation of magnetic 
monopole pairs, which would drive the system from our 
Coulomb phase to a phase in which the U(l) gauge field is 
gapped. Since the energy of a Hopfion is essentially mag- 



netic field energy according to the discussion in Sections [HD 
and [IVJ monopole pair condensation should leave the Hop- 
fion gapless and, perhaps, even soften its dispersion relation 
by suppressing magnetic field energy. Therefore, this does 
not lead to a fully-gapped phase and, as a result of the soft- 
ened spectrum, may even lift the zero modes. Within the con- 
text of the fermion dimer model, the natural approach is to 
consider a non-bipartite lattice. After all, in two dimensions, 
the dimer model on a triangular lattice has a gapped liquid 
phase^^. A possible three dimensional lattice to consider is 
the stacked triangular lattice. This lattice can be chosen to 
have TT-flux in each plaquette as required. Although there may 
be problems with this approach (for example, we could find 
that the monomers become confined in planes, or that there 
is no liquid phase), it may allow the construction of a gapped 
three dimensional model and quantum Monte Carlo simula- 
tions of this system would be useful. We have not been able 
to construct a field theoretic version of this model, however. In 
this model, suppose the Hopfion is a gapped, localized object. 
In this case, implementing an operation such as twisting both 
quasi-particles would require twisting each quasi-particle sep- 
arately, emitting a Hopfion pair, and then bringing the Hopfion 
pair together and annihilating them. However, it seems unnat- 
ural for the Hopfion to be gapped: suppose a quasi-particle 
is dragged along a long closed loop and returns to its origin 
with a changed sign, causing a Hopfion to be emitted. At 
what point along this path does the system have to pay an en- 
ergy price to emit the Hopfion? Since the system is gapped 
and liquid-like, all points on the path should be the same, so 
why should one point be distinct? One possible resolution is 
that the Hopfion is gapped in the bulk, but that there is a zero 
energy bound state of a Hopfion and a defect, reminiscent of 
what happens in the triangular lattice model in two dimen- 
sions. Again, this is a problem for quantum Monte Carlo to 
study. 
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Appendix A: Relation between the Hopfion and the gauge field 
instanton in the (4+l)-dimensional theory 



In this appendix, we will present a more detailed discussion 
of the relation between the Hopfion of the order parameter and 
the gauge field instanton in the (4+l)-d theory discussed in 



Sec. IV B In the interest of generality, we consider the generic 
action in Eq. ( [73] ) with N copies of Weyl fermions forming the 
fundamental representation of SU(N), coupled to an SU(N) 
gauge field and a Higgs field n''^ e SU{N)/SO{N). 

We first consider the symmetry unbroken phase n^^ = 0. 
The axial anomaly leads to the equation 



5^i^ 



1 



327r2 



Tr [e'^^'^^U^Ur] 



with /^j, the SU(N) gauge curvature, and 



f = c^a^c 



(Al) 



(A2) 



is the number current of the Weyl fermion. In the 4 + 1-d 
regularization of this theory, the Weyl fermions live on the 



boundary. For the slab geometry shown in Fig. (IV B), an 
instanton configuration of f^jy with the second Chern number 



Id' 



^3271 



,Tr [e^^^^/^^/^^] = TV pumps TV fermions 



C2 

from one boundary to the other. The instanton configuration 
is determined by the homotopy group 7Ts{SU{N)) = Z. Con- 
sider an instanton configuration with C2 = N in the (bound- 
ary) space-time R^ with the boundary condition /^j, -^ for 
time t -^ ±00. If we choose a gauge so that the gauge po- 
tential a^ is continuous and a^{-K^t) -^ for t -^ —00, 
then a^{-K^t) -^ —ig~^{yi)d^g{yi) for t -^ +00, in which 
^(x) G SU(N) defines a mapping x -^ ^(x) with nontrivial 
winding number 

^=^1 ^'xe^^'^Tr [g-'d,gg-'djgg-'dkg] (A3) 



Now we consider the Higgs phase with the order parame- 
ter n^^{-K^t) defining a map from the space-time to the tar- 
get space SU{N)/SO{N). A Hopfion is a soliton configura- 
tion n^^(x, 0) with nontrivial winding number determined by 

TTs {SU{N)/SO{N)) = Z2. For unitary symmetric matrices 



-,ab 



satisfying n 



1, one can always find a gauge 



transformation ^(x) G SU(N) such that 
n(x) = g^g{x.) 



(A4) 



The choice of ^(x) is not unique since for any 0(x) G 
SO(N), ^(x) = 0(x)^(x) also satisfies n(x) = g^g{x.). 
Thus ^(x) is determined modulo a SO{N) factor multiplied 
from the left, which is consistent with the fact that the or- 
der parameter space is the coset SU(N)/SO(N). The Z2 
winding number of the configuration n(x) is related to the 
integer- valued winding number of the gauge transformation 
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^(x) given in Eq. ( |A3[ ) by the following equation: 

Nz^ [n(x)] = N [^(x)] mod 2 



For n = g^ g one obtains 



-ig ^d^g 



(All) 
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mod 2 



(A5) 



As is discussed above, ^(x) has an SO(N) ambiguity. For 
two configurations ^(x) and g(x) = 0(x)^(x) with 0(x) G 
SO(N), the winding number ( A3 ) is different by 



7V[^(x)]-7V[^(x)] = A^[0(x)] 



Thus we see that a Hopfion configuration with nonzero Nz^ 
defined in Eq. ( |A5| ) is always associated with a pure gauge 
configuration with odd winding number N [^(x)]. The Hop- 
fion number is invariant under smooth deformation of the field 
n(x), and thus remains conserved in time evolution unless the 
orderparameter vanishes. If the order parameter is allowed to 
vanish, such as by pair creation of monopole pairs, the Hop- 
fion number can be changed. An example of such process 



-Q-^djOQ-^dkO] 



(A6) 



N [0(x)] is the winding number of the orthogonal matrix 
0(x), which is always an even number. Thus the winding 
number Nz^ [^(x)] = ^ [^(^)] "^od 2 is not affected by the 
choice of ^(x) and is determined completely by n(x) = g^ g. 

The relation between the winding numbers of n(x) and 
^(x) presented above can be verified directly, but we also 
would like to provide an alternative explanation in a more 
mathematically rigorous way. The facts we discussed here 
arise from the exact sequence 

1 -> SO{N) -^ SU{N) -^ SU{N)/SO{N) -^ 1 (A7) 

which leads to the long exact sequence 



is illustrated in Fig. (17). Consider a generic Hopfion cre- 
ation process, which connects the initial configuration rii (x) 
with Hopfion number Nz2 = and the final configuration 
712 (x) with Hopfion number Nz2 = 1- Independent from the 
detail of the Hopfion creation process, according to the dis- 
cussion above one must simultaneously create an instanton of 
the gauge field. Defining rii = ^f ^i(x), n2 = ^2^^2(x), one 
obtains 



C2 = 



1 
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Tr [e^^^V^./ar] = N [g,{^)] - N [^i(x) 



(A12) 



Consequently 

C2 mod 2 = Nz, [n2(x)] - Nz, [ni(x)] = 1 (A13) 

Thus we see that the gauge field configuration associated with 
this Hopfion creation process is an instanton with odd Chern 
number, which corresponds to a change of the Weyl fermion 
number 



^ ^n+i {SU{N)/SO{N)) ^ TTn {SO{N)) ^ TTr, {SU{N)) 

^ TTn {SU{N)/SO{N)) ^ TTn-l {SO{N)) ^ ... (A8) 



ANf = C2 ^ ANf mod 2 = 1 



(A14) 



Considering that 7r2{SO{N)) = 1 and 

7r4{SU{N)/SO{N)) = 1 (with 1 standing for a trivial 
homotopy group), we have 



ns {SO{N)) ^ ^3 {SU{N)) 
TTs {SU{N)/SO{N)) -^ 1 



(A9) 



From this exact sequence we see that the homotopy 
group TTs {SU{N)/SO{N)) = Z2 is a quotient group of 
TTs {SU{N)) = Z and tts {S0{N)) = Z. In other words there 
is an injective projection from tts {S0{N)) to tts {SU{N)) 
which maps each SO(N) configuration with winding num- 
ber n G ZtoaSU(N) configuration with winding number 
2n. The nontrivial configurations of the coset SU(N)/SO(N) 
are the image of the SU{N) configurations with odd winding 
number. 



With the preparation above, now we consider the gauge 
field coupling to the Higgs field. With the configuration n(x) 
of the Higgs field, the gauge field configuration in the ground 
state is a pure gauge determined by the condition 



D^n 



d^n 



■ T 

ia,,n 



ma,. 







(AlO) 



Thus the fermion number of the boundary system changes by 
an odd number. Although in the Higgs phase the Fermion 
number conservation is explicitly broken, the Fermion number 
modulo 2 is still conserved in the classical level. 

In summary, the discussion above leads to the conclusion 
that the Z2 symmetry of Fermion number modular 2 is bro- 
ken anomalously for this boundary system, and the instanton 
configuration that produces this anomaly is a Hopfion creation 
process. One can also provide a complementary picture in the 
bulk of the (4 + l)-d regularized theory. To understand the 
relation between Hopfion number and fermion number par- 
ity in the bulk, we introduce a /7(1) probe gauge field Aj 
(I = 0, 1, .., 4) and couple the fermions to the field Aj + a/. 
Integrating over the fermions leads to a Chern- Simons term of 
the gauge field Aj + a/ (c.f. RefJ5j: 

5eff = ^Jd^xe'^'^'^^TTiiAj^aj) 

■dj {Ak + aK) Ol {Am + au) + ...] (A15) 

with . . . standing for the terms with higher orders of Aj + a/ 
required by gauge invariance. The fermion number current 
couples to the gauge field Ai and is described by the response 
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equation 



SS, 



eff 
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JJKLM 



A,=0 
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Tr [/jk/lm] (A16) 



In any 4-d spatial region D, the fermion number is 



NpiD) = J fd'x = I d'x^e'^^^^Tv [/jk/lm] 



(A17) 



In the Higgs phase, the fermion number conservation is bro- 
ken but the fermion number mod 2 is still conserved. Thus the 
equation above still holds modular 2. 

In the 4-d space the Hopfion is a point-like defect. The 
process of creating of Hopfion on the boundary can also be 
thought as pulling a Hopfion defect across the surface into the 



bulk, as is illustrated in Fig. (17 ). Consider a spatial region 



D including a Hopfion defect and its neighborhood. As long 
as the region D is large enough to cover the core region of the 
Hopfion, the gauge field at the boundary S = dD is a pure 
gauge. By definition, the configuration n(x) on E has an odd 
winding number Nz2 [^(x)] = 1, which corresponds to a pure 
gauge configuration with an odd winding number Nj^ defined 



in Eq. ( A3 ). Thus the fermion number in Region D satisfies 



N, 



<D) = f 
Jd 



<fx- 



1 



^OJKLM 
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Tr [/jkIlm] = N^ mod 2 



(A18) 



which demonstrates that each Hopfion carries an odd fermion 
number. When the Hopfion is created and then moved across 
the bulk to the other surface, one fermion (or odd number of 
fermions) is pumped from one boundary to the other. 
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